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Abstract 

We study analytically, via the Newman-Penrose formalism, the late time 
decay of scalar, electromagnetic, and gravitational perturbations outside a 
realistic rotating (Kerr) black hole. We find a power-law decay at timelike 
infinity, as well as at null infinity and along the event horizon (EH). For 
generic initial data we derive the power-law indices for all radiating modes 
of the various fields. We also give an exact analytic expression (accurate 
to leading order in \/t) for the r dependence of the late time tail at any 
r. Some of our main conclusions are the following, (i) For generic initial 
data, the late time behavior of the fields is dominated by the mode / = |s| 
(with s being the spin parameter), which dies off at fixed r as t"^!'*!"^ — 
as in the Schwarzschild background, (ii) However, other modes admit decay 
rates slower than in the Schwarzschild case, (iii) For s > fields, non-axially 
symmetric modes dominate the late time behavior along the EH. These modes 
oscillate along the null generators of the EH. 
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I. INTRODUCTION 



Small perturbations of the Schwarzschild black hole (SBH) geometry die off at late time 
with an inverse power-law tail. This well-known phenomenon was discovered by Price early 
in the 1970s. Price explored the dynamics of linear scalar and metric perturbations 
(and that of all integer-spin fields in the Newman- Penrose formalism p[) propagating on the 
SBH background. His analysis provided a detailed description of the relaxation mechanism 
through which the black hole (BH) exterior settles down at late time into its stationary "no 
hair" state. In particular, Price was able to characterize the actual form of the late time 
fall-off of the perturbations: He found that any radiative multipole mode /, m of an initially 
compact linear perturbation dies off at late time as (where t is the Schwarzschild time 

coordinate). In the case there exists an initially static multipole mode l,m it will decay as 
^-21-2^ These power-law decay tails were found to be the same for all kinds of perturbations, 
whether scalar, electromagnetic or gravitational (and in this respect, the scalar field model 
proved to be a useful toy model for more realistic fields). 

Price's results were later reproduced using several different approaches, both analytical 
and numerical and were generalized to other spherically symmetric BH spacetimes 

PJ9|-[T^. (A brief review of the works on this subject can be found in the Introduction of 
Ref. [0.) The validity of the perturbative (linear) approach was supported by numerical 
analyses of the fully nonlinear dynamics [p!5| , pT| , indicating virtually the same power-law 
indices for the late time decay. 

For a scalar field on the background of a SBH, power-law decay tails were found to be 
exhibited also at future null infinity and along the future event horizon @,||. It was 

shown that at null infinity the scalar field dies off with respect to retarded time u as 
(for a compact initial mode) or as (for a static initial mode). The decay of the scalar 

perturbation along the event horizon (EH) was found to be ti'^'-s qj, -y-^'-s ^f^j, ^ compact 
or a static initial mode, respectively), where v is the (Eddington-Finkelstein) advanced-time 
coordinate. 

As was already explained by Price, the late time tails of decay outside spherically sym- 
metric BHs originate from backscattering of the outgoing radiation off spacetime curvature 
at very large distances. In the framework of a frequency- domain perturbation analysis 
PJ5IJTII , these tails are explained in terms of a branch cut in the (frequency domain) Green's 
function, which is, again, associated with the form of the curvature-induced potential at 
large distance. This suggests that the form of the decay at late time reflects (and is affected 
by) merely the large distance structure of spacetime, and may be independent of the exis- 
tence or absence of an event horizon. This assertion found further support in the analysis 
by Gundlach et al. [O] , who studied the purely spherical collapse of a self-gravitating mini- 
mally coupled scalar field. It was demonstrated numerically that in this case late time tails 
form even when the collapse fails to create a black hole. (On the other hand, quasi-normal 
ringing are found to dominate the early stage of the waves' evolution only if a BH forms.) 

Until quite recently, the issue of the late time decay of BH perturbations has been consid- 
ered only in spherically symmetric models of BHs. It is known, however, that astrophysically 
realistic BHs are spinning , and thus are not spherically symmetric but are rather of the 
axially symmetric Kerr type. Moreover, it is suggested, in virtue of the "no hair" principle 
p!7| , that the Kerr black hole (KBH) might be the only realistic BH (realistic BHs are not 
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expected to carry a significant amount of net electric charge). Hence, generalization of the 
above mentioned analyses to the KBH case seems to be of an obvious importance. Still, 
such a generalization has awaited almost three decades, as the lack of spherical symmetry 
in the Kerr background makes both analytical and numerical exploration significantly more 
complicated. 

The "no hair" principle for BHs implies that perturbations of the KBH must "radiate 
away" at late time. No further information is available from this general principle as to the 
actual details of this decay process. A question arises as to what effect rotation has on the 
form of the late time tails. More basically, does the decay of the perturbing fields still obey 
a power law? If so, are the power indices the same as in the SBH case? 

Such questions were addressed only recently, by several authors. First, Krivan et al. 
carried out a numerical simulation of the evolution of linear scalar and gravitational p9[] 
waves on the background of a Kerr black hole. The first analytic treatment of this problem, 
for a scalar field, was later presented by Barack and Ori [^,^ (following a preliminary 
analysis by Ori [Q). Most recently. Hod used a different approach to study the late time 
decay of both a scalar field and nonzero-spin Newman- Penrose fields p^ , |25[] on the Kerr 
background. (Hod's analysis, which follows some preliminary considerations by Andersson 
1^, is carried out in the frequency domain, whereas Barack and Ori use a time domain 
analysis.) Finally, the analytic progress has motivated a further numerical study by Krivan 

mm- 

The above analyses all indicate that power-law tails of decay are exhibited in the Kerr 
background as well. In the Kerr case, however, the lack of spherical symmetry couples 
between various multipole modes, which results in the power-law indices of specific modes 
being found to be different, in general, from the ones obtained in spherically symmetric BHs. 
Another interesting phenomenon caused by rotation (first observed in l^^j) is the oscillatory 
nature of the late time tails along the null generators of the EH of the Kerr BH for nonaxially 
symmetric perturbation modes. 

The purpose of the present paper is to extend the analysis described in Ref. |21[] to 



electromagnetic and gravitational perturbations of the Kerr background, and supply the full 
technical details of our approach. (In Ref. |^ we merely outlined the application of our 
technical scheme to a scalar field, and gave a brief description of the results in this case.) 
The analysis to be described in this paper provides a more complete and accurate picture 
of the late time decay of physical fields, than already available. Among the results which 
appear here for the first time: 

(i) We derive the form of the late time tail for all radiative modes anywhere outside the KBH 
(i.e. at all distances). We also give an exact analytic expression for the radial dependence 
of this tail. (In p3|-|25| Hod only analyzes the asymptotic behavior at very large distance, 
and along the EH.) 

(ii) A careful analysis of the decay along the EH reveals an interesting phenomenon: For 
s > fields, it is the oscillatory nonaxially symmetric (m > 0) modes which dominate 
the late time behavior there. This result has important implications to the structure of 
the singularity at the inner horizon of the KBH p9[. In a different paper p8| Barack and 
Ori further explore and explain this phenomenon, and discuss the reason for the incorrect 
prediction made in Ref. for the decay rate along the EH in the s > case. (Recently 

following the appearance of Ref. |^ as report no. ^r-qc/9907085| . Hod has corrected 
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his result.) 

For simplicity, we shall refer in this paper only to the (most realistic) situation in which 
the initial perturbation has a rather generic angular distribution, such that it is composed 
of all multipole modes and, in particular, the lowest radiatable one. It will be explained why 
an initial setup in which the lowest modes are missing is more complicated to explore using 
our approach. Hod's analysis |]23|-p5[] provides predictions for this well; however, 

these are not in agreement with recent numerical results by Krivan It thus seems that 
further analytic work is needed (using either Hod's method or ours) to clarify this point. 

The arrangement of this paper is as follows. In Sec. ^ we briefly review the subject of 
perturbations of the Kerr geometry via the Newman-Penrose formalism. We then reduce 
the master perturbation equation in the time domain to obtain a coupled set of time-radial 
equations for the various modes. The evolution problem for each of the modes is mathe- 
matically formulated as a characteristic initial value problem. In Sec. |T| we analyze the 
late time behavior of the various modes at null infinity. To that end we apply the iterative 
scheme, which is, basically, an extension of a technique previously tested for a scalar field in 
the SBH background . In Sec. ^ we introduce the late time expansion (LTE) scheme, 
which allows a global treatment of the decay at late time. The late time behavior of fields at 
a fixed distance outside the EH is dominated by the leading term of the late time expansion, 
for which we derive an analytic expression in Sec. 0. In Sec. |VI| we then carefully explore 
the behavior along the EH itself. Finally, in Sec. |VII| , we use the LTE scheme combined 
with the results at null infinity, in order to derive the late time decay rates for all modes of 
the fields at any fixed distance. We conclude (in Sec. [VI I H ) by summarizing our results and 
discussing their relation to other works. 

II. MODE-COUPLED FIELD EQUATION 

A. Perturbations of the Kerr geometry via the Newman-Penrose formalism 

(definitions and a brief review) 

The line element in Kerr spacetime reads, in Boyer-Lindquist (BL) coordinates t, r, 6, if, 

ds^ = - (1 - 2Mr/S) dt^ + (S/A) dr^ + JldO^ + (r^ + + 2a^Mr sin^ O/T.) sin^ 6 dip^ 
- (AaMr sin^ ^/s) d'p dt, (1) 

where M and a are, correspondingly, the BH's mass and specific angular momentum, E = 
+ cos^ 6, and 

A = - 2Mr + a^. (2) 

(Throughout this paper we use relativistic units, with c = G = 1.) We shall consider in this 
paper only a BH solution with \a\ < M: the extremal case, \a\ = M, requires a separate 
treatment, as we later briefly explain. The event and inner horizons of the (non-extremal) 
KBH are the two null 3-surfaces r = r_|_ and r = r_, respectively, where 

r± = M ± VM2 - a2 (3) 
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are the two roots of the "horizons function" A(r) defined in Eq. (^. 

To discuss perturbations of the Kerr spacetime via the Newman- Penrose (NP) formahsm 
||3CI|| , introduce Kinnersley's null tetrad basis |31[] [l'^, n^, m^, m*^] (where an asterisk denotes 
complex conjugation). In BL coordinates, the "legs" of this tetrad are given by 



P = A" 



r^ + a^A,0, a 



(2S)~i p + a2,-A,0,a 
(2^/2p)-i [iasin^,0,l,Vsine] 



(4) 



(with the fourth tetrad vector obtained from by complex conjugation), where p = 
r + ia cos 6. In the framework of the NP formalism, the gravitational field in vacuum is 
completely described by five complex scalars \l/o, . . . , ^1/4, constructed from the Weyl ten- 
sor Cajj-yS by projecting it on the above tetrad basis. Likewise, the electromagnetic field is 
completely characterized by the three complex scalars ipQ,ipi,ip2, constructed by similarly 
projecting the Maxwell tensor F^,^. In particular, 

^0 = -Ca^^srm'^Pm^ and ^4 = -C^f^^sn"m*^n^m*^ (5) 

represent, respectively, the ingoing and outgoing radiative parts of the Weyl tensor, and 



and 



^2 



(6) 



represent the ingoing and outgoing radiative parts of the electromagnetic field. 

In the (unperturbed) Kerr background all Weyl scalars but \l/2 vanish (as directly im- 
plied by the Goldberg-Sachs theorem in view of Kerr spacetime being of Petrov type D; 
see Sees. 9b, 9c in [0). In the framework of a linear perturbation analysis, the symbols 
^'o, ^1, ^^2, ^3, ^4 and ipo,6ipi,ip2 are thus used to represent first-order perturbations of 
the corresponding fields (with 6"^ 2 = "^2 — vj/^'^'^'^s™"'^'^^ etc.). One can show (see Sec. 29b 



in [^) that \l/o and \l/4, and also v^o and ip2, are invariant under gauge transformations 
(namely, under infinitesimal rotations of the null basis and infinitesimal coordinate trans- 
formations). The scalars and \I'3 are not gauge invariant, and may be nullified by a 
suitable rotation of the null frame. The entities 5\l/2 and 6(pi represent perturbations of 
the "Coulomb-like", non-radiative, part of the fields (in fact, one can also nullify 5\E'2 by a 
suitable infinitesimal coordinate transformation.) It is therefore only the scalars defined in 
Eqs. (^ and (|^) which carry significant information about the radiative part of the fields. 
(Note, however, that gauge invariance of the radiative fields is guaranteed only within the 
framework of linear perturbation theory.) 

Teukolsky first obtained a single master perturbation equation governing linear per- 
turbations of scalar, electromagnetic, and gravitational fields. In vacuum, this master per- 
turbation equation reads 



A 



sin^ 9 



4Mar f LJ\ iu^ 

+ ^ '^.*^+Ia sin^J 



sin 9 

+ (s2 cot^^-s)^ 



sin^) „ - 2s 
0. 



a{r — M) ^ i cos 9 



A 



sin'^ 



M(r2 



A 



— r — ia cos 9 



(7) 
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Here, r, 5, y?) represents the various radiative fields according to the following list: 

<|. = ^^=o (scalar field), 

^4 = (p*)-4^^=-2. (8) 

The master equation (0) is fully separable only in the frequency domain, by means of the 
spin-weighted spheroidal harmonic functions S^^"^{—a'^u^, cos 6) ||33|, where u is the temporal 
frequency (the separated equations are referred to as "Teukolsky's equation"). Because the 
functions S"''™(^) are u dependent, separation of the 9 dependence is not possible in the 
time domain (namely, without first decomposing the field into its Fourier components). 



B. Reduction of the master field equation in the time domain 

The target of the present work is to explore the behavior of the fields at late time. 



In principle, the analysis can be carried out in the frequency domain, as in Refs. |23H25 
In this technique, the requested temporal behavior is finally to be extracted by an inverse 
Fourier transform. Our analysis is based on a different approach, motivated by the following 
argument: In the late time, stationary, limit {t oo), one expects the very low frequency 
(a; — > 0) Fourier modes to dominate the behavior. For such waves, the functions e«™-¥'5'«'™(^) 
reduce to the spin-weighted spherical harmonic functions ¥'^''"^{9, ip) This may motivate 
one to try and extract the angular dependence of the fields by using the functions y'"*^. 
As a result of the lack of spherical symmetry, the resulting (time-domain) field equations 
will possess coupling between the various multipole modes I; however, one should expect 
this coupling to be "small" , in a sense, at late time. In the sequel we show how this coupling 
can be treated in an iterative manner, in both the frameworks of the iterative expansion 
(Sec. pT[) and the late time expansion (Sec. [TVl) . 

Led by the above argument, we expand the fields as 

oo / 

*^(t, r, 9, ip) = (r^ + a')-i/2A-^/2 ^ ^ ¥''""{9, ip) ^^''"(t, r), (9) 

;=|s| m=—l 

where the radial factor in front of the summation symbols is introduced for convenience 
[as it eliminates the term oc ip^y, from Eq. (p3D below]. Note that the summation over 
modes / excludes the / < |s| modes, which are nonradiative (for a discussion regarding the 
nonradiatable modes, see Sec. HID 4 in Ref. 0). 

Inserting the expansion (^ into the master field equation (|^), we obtain 

^ y^''"(^, if) [b{t, r)V"'" - sin^ 9 (V'''"),^ + 2ms cos 9 (V-''™) J = 0, (10) 

Im 



where D(t,r) is a certain differential operator independent of 9, (p. Note in Eq. (ITU]) how 



the two last terms in the squared brackets (proportional to a) avoid a full separation of 
variables. 
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Now, the product cos 9 ■ y*'"^ can be re-expanded in terms of the functions y (which 
form a complete set of functions on the unit 2-sphere for each s). The "matrix elements" of 
the function cos 9 with respect to the y basis are given by [^] 

r / 2/ + 1 \ 

{si' m\ cos 9\slm) = j dn iY'^"")* cos 9 {Y'^"") = l ^^l' + i ] {^'^rnO\l'm){ll - sO\l' - s), 

(11) 

where (jij2^i^2|j^) are the standard Clebsch-Gordan coefficients [^. We find that 

cos^ ■ y' = c'+'r'+' + 4y' + c'-^y'~\ (12) 

where 



i2(2«-l)(2/+l) 
ms 
1(1+1)' 



(13) 



(Here, as we shall often do below, we omit the indices s,m for the sake of brevity.) This 
also easily leads to 

- sin^ 9-Y' = cLt' y'+2 + ci+i y'+i + + c'-' y'-i + c'-^ y'-\ (m) 

where 

_ J+iJ 

= <(4+i + c[,), 

= (cL)' + {c\f + (4)2 - 1, (15) 



(ji — J-iJ 



All constant coefficients c' and C' are nonvanishing, with only the following exceptions: 

(i) Ci and cL vanish for Z = |m| or / = 

(ii) C\ vanishes for |m| < / < |m| + 1 or |s| < / < |s| + 1. 

(iii) Cq, C^, and CL vanish for m = or s = 0. 

(It can be verified that C\ is always negative definite.) 

Substituting Eqs. (0) and (|1^) in Eq. (^) we obtain by the orthogonality of the functions 

"yslva 

i)(t,r)V''+X(^/^'±\V'''^2) ^ (16) 

for each l,m,s satisfying / > \m\ and I > \s\, where D is yet another differential operator, 
and X is a functional describing coupling between the I mode and the I ± 1 and / ± 2 modes: 

T(V^'±\ ^'±2) ^ ^2 (^^^^^^+2 + + cl^i-i + + 2ms (c^^'+i + c'_V^'-i) ^ . 

(17) 
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Note that, obviously, modes of different m do not interact with each other, as the Kerr 
geometry is axially symmetric. We also point out that the scalar field case is special, in that 
for this case the / mode does not interact with the / ± 1 modes, but only with the / ± 2 
modes (recall that the interaction coefficients vanish for s = 0). 

To write Eq. (|1^) explicitly in a convenient form, we introduce the advanced and retarded 
time coordinates, defined, respectively, by 



V = t + and u = t — 



(18) 



(which are nonetheless not null coordinates in Kerr spacetime). Here, the "tortoise" radial 
coordinate r* is defined by 



r^. = r — r+ + (2k+) ^ln2;+ — (2k_) ^lnz_ 

with K,± being the horizons' "surface gravity" parameters, 

r_|_ — r_ 

and where the dimensionless radial variables z± are given by 

r — r± 



r I — r_ 



(19) 



(20) 



(21) 



[Note the relation A = z+z_(r_|_ — r_)^. Also, recall that we are dealing in this paper only 
with non-extremal black holes, for which r+ > and, consequently, z± are well defined.] 
The coordinate r^,, satisfying dr^,/dr = (a^ +r^)/A, increases monotonically with r from 
— oo (at the EH) to +oo (at spacelike infinity). Later we shall find useful the asymptotic 
relations 



(r-i-— r_)2 



2+ (for <C — M, r r+), 
(for ^ M, r — > oo). 



The explicit form of the mode-coupled field equation ([T6|) is nowQ 
where V"'(r), i?'(r), and K{r) are radial functions given by 



(22) 



(23) 



AV^(r) = (r^ + a^)-^ 



{I- s){l + s + 1)A - m^a-^ - 2isma{r - M) 



d r 



dr dr 



(24a) 



2RUr) = (r2 + a2)-2 



2imMar - sM{r'^ - a^) + sA(r + iacg) 



(24b) 



^ In Eq. (23), the derivative du is taken with fixed v, dy is taken with fixed n, and dt with fixed r. 
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4:K(r) 



(r' + a^)-^A. 



(24c) 



Of importance will be the large-r asymptotic forms of these functions, which in terms of 
the expansion read 



4r? 



2r2 



+ 



(Inr*)" 



(25a) 



lacr 



-3M + 2M ln(- 



2r2 



+ 



(Inr^ 



(25b) 



ir(r) 



1 

4^2 



M 



1 + 2 In ( 



2r3 



(Inr*) 



(25c) 



Equation (23) already provides a qualitative picture of one aspect of the fields' evolution: 
Unlike in spherically symmetric spacetimes, multipole modes of different / interact with 
each other while propagating on the Kerr background.^ For example, if a physical (s 7^ 0) 
perturbation is initially composed of a pure l,m mode, then, in general, we may expect all 
possible modes l',m — namely, all modes with /' > max(|m|, |s|) — to be generated while the 
perturbation evolves in time. [The unrealistic case of a scalar field (s = 0) is special, as for 
this field only those of the above modes with even / — I' are expected to be excited.] Note 
that all interaction terms in Eq. (|23|) involve derivatives of the field with respect to t. It is 
this feature which, by means of the LTE scheme, allows one to effectively decouple Eq. (p3[), 
as we show in Sec. iT 



C. The initial setup 

We shall consider an initial perturbation in the form of a compact outgoing pulse of 
radiation, which is relatively short, yet arbitrarily shaped. We take this pulse to be emitted 
a.t u = Uq, V = Vq, and without limiting the generality we take, for simplicity, Vq = 0. 
We further assume that the initial pulse has a rather generic angular shape, so that it is 
composed of all multipole modes l,m [and in particular, for each m it contains the lowest 
possible mode, / = max(|m|, |s|)]. Formulated mathematically, this initial setup takes the 
form0 



2 See, however, our remark in the concluding section, with regard to the definition of "multipole 
moments" in the Kerr geometry being somewhat ambiguous. 

^ It should be noted that, strictly speaking, these initial data [supplemented to the coupled 
field equation (23)] do not form a well posed characteristic initial- value problem, as u =const 



and V =const are not characteristic hypersurfaces of the Kerr geometry (these hypersurfaces are 
timelike rather than null). We further comment on this issue below. 
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(26) 



where, for each / and m, r''"(-u) is an arbitrary (but nonvanishing) function with a compact 
support between u = Uq and (say) u = Ui > Uq, with Ui — uq \uo\. This type of initial 
data corresponds to the physical scenario in which no ingoing radiation is coming from past 
null infinity. 

It will be assumed in the following that the initial pulse is emitted far away from the BH; 
namely, we take —uq ^ M. This assumption greatly simplifies our analysis (as we explain 
in Sec. |Tl|; cf. 0,§]); yet, it seems reasonable to expect the late time behavior in this case 
to remain characteristic of the general situation. 



III. LATE TIME BEHAVIOR AT NULL INFINITY: THE ITERATIVE SCHEME 

In this section we derive the form of the late time decay at future null infinity, that is 
for f — > oo at finite u ^ M. This has two main motivations: First, the results at null 
infinity will appear to serve in the framework of the late time expansion scheme as necessary 
"boundary conditions" for the global late time evolution problem (as we discuss in Sec. |rv| ). 
The results at null infinity also have their own physical significance, for the following reason: 
Consider a static (fixed-r) observer located at very large distance. Let Am and At represent, 
respectively, the retarded and the static observer's time elapsed since this observer gets the 
first signal from the perturbing field. With respect to this observer, the relevant information 
about the decay at finite Am = At is the one calculated at the null infinity domain (which 
we may call, in this context, the "astrophysical zone" of the waves). Only when the time 
lapse becomes infinitely large (while r remains finite) does this observer enters the "future 
timelike infinity" zone, t ^ r (the late time behavior in this domain will be discussed in 
Sec. IVITD. 



For this part of the analysis (namely, for the derivation of the late time tails at null 
infinity) we apply the iterative scheme, first developed and tested for scalar waves on the 
Schwarzschild background in Refs. Since the technical details of the calculations in- 

volved are often very similar to those in the above references, we mainly describe here the 
new results (for spin-s fields in Kerr), and direct the reader to Refs. for further details. 

A. Formulation of the iterative scheme 

We define 

fio(r) ^ ^ (27) 



and 



SR\r) = R\r) - Ro{r). (28) 



10 



The functions V^'(t') and Ro{r) extract at large r the (flat space) asymptotic behavior of the 
functions V^'(r) and i?'(r) appearing in Eq. (^31) . The "curvature induced" residual part of 
these functions is represented by 6V'' (r) and SR'^r). 

We now decompose each of the functions -0*'™ (for each l,m) as 

^l;' = 4 + ^[ + ^l + ..., (29) 
such that each of the functions ipl^ satisfies the field equation 

^i,™ + V^ir)^pl + RoM, = (30) 

with S'^^Q = and with 

Sn>o = - SR\^Pl.,),t - Kir) [a'Cl, i^l_,),u + ^i-'i)] • (31) 

We take the initial conditions for the various functions ipl^ to be V'i=o — V"' and ip^^Q = on 
the initial surfaces u = uq and f = 0. 

Equation (POD, supplemented by the above initial data, constitutes a hierarchy of charac- 
teristic initial data problems for the various functions ip'^. Formal summation over n recovers 
the original evolution problem for ipK For each n > (and for all /), the function ip'^ admits 
an inhomogeneous field equation, with a source term depending only on the functions V'n'<n 
preceding it in the hierarchy, and with the function V'i=o satisfying a closed homogeneous 
equation. This structure allows one, in principle, to solve for all functions ip'^ in an iterative 
manner: first for all modes / of 'i/'i=0' which then serve as sources to '4'^n=ii l^ general, 
each function ip'^ shall have sources coming from the modes Z, / ± 1, and / ± 2 of ipn~i- 

Of course, the effectiveness of the proposed iteration scheme crucially depends on its con- 
vergence properties. In that respect, the scheme as formulated above may seem problematic, 
because, while the "zeroth order" (n = 0) field equation well approximates the actual field 
equation (^) at large distance, it fails to do so at the highly curved small-r region (actually, 
the functions Vq and i?o; as defined above, diverge at = 0). 

It is possible (by redefining Vq and Rq at small r) to construct a more sophisticated 
iteration scheme that would account for the small-r region of spacetime as well — as was 
done in Ref. for the case of a scalar field in the Schwarzschild spacetime. In that case, 
it was demonstrated that the small-r details of the background geometry have merely 
a negligible effect on the late time behavior at null infinity. In the Schwarzschild case, we 
were able to greatly simplify the analysis by considering a toy model of a thin spherically 
symmetric shell of matter having a fiat interior. In that model, the late time behavior of 
the scalar field at null infinity turned out to well approximate the actual behavior in a 
"complete" Schwarzschild model (see [|^|[ for details). 

The above results all indicate that the late time decay of the scalar field is predominantly 
governed by the large-r structure of the Schwarzschild spacetime. (This conclusion stemmed 
already from several previous works, as we mentioned in the Introduction.) We shall assume 
in this paper that the same is also valid for all fields propagating on the Kerr background. 
To be concrete, we will consider in this section, only for the sake of the calculation at null 
infinity, a model in which the Kerr interior geometry is replaced by a fiat (Minkowski) 
manifold. There is no way of smoothly attaching a fiat interior to a Kerr exterior through 
a thin spherically symmetric material shell (as was done in in the Schwarzschild case); 
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alternatively, we take this attachment to be made through a material layer of a finite width. 
The external "radius" of this layer should be of order of a few M. Based on our experience 
with the Schwarzschild case, we expect the details of our model at small distances (and, in 
particular, the internal structure of the material layer) to merely have a negligible affect on 
the behavior of the waves at null infinity, at u ^ M . 

We emphasize that the above simplified model is adopted only for the analysis at null 
infinity, where it greatly reduces the amount of technical details one should deal with. For 
the global analysis carried out in the rest of this paper, the "complete" KBH geometry shall 
be considered. 



B. Zeroth order iteration term 

By definition, the function ipl^^Q obeys (for all /) the homogeneous equation 

C + K,'(r)^^ + i?o(r)<, = 0, (32) 

with the initial conditions iP\^^q{v = 0) = V\u) and ipl^^^iu = uq) = 0. Equation (|32|) simply 
describes the free propagation of the field ipn=o Minkowski spacetime, provided that is 
replaced with the radial Minkowski coordinate r. 
The general solution to this equation readsQ 

/general _ ^ ^sl 9o\u) , ^~sl ho\v) , , 

~ho 'j^^^^uj^^h ' (^7^^' ^^^^ 

in which QqIu) and hQ{v) are arbitrary functions (with their parenthetical superindices indi- 
cating the number of differentiations), and where A''^ are constant coefficients given by 



= Jr /\.r (34) 



The above initial data for il'L=o uniquely determine a specific solution for this function [but 
not for either of the functions go{u) and Hq^v) in separate). For the outgoing pulse initial 
setup it is possible to write this specific solution in a convenient form by taking hQ{v) = 0, 
in which case we have 

i'L.-pfj^. (35) 
with the function goiu) given byQ 



Note that the homogeneous equation (^) is invariant upon simultaneously transforming u ^ v 
and s ^ —s. The two sums in Eq. ( p^ ) constitute two independent homogeneous solutions, which 
are obtained from each other by this transformation. 

^ The calculation leading to Eq. (|3^) is a straightforward generalization of the one described in 
detail in Ref. for s = 0. 
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g'>\u) = e{-u) X r{u/uT'^\u - uT'-\-urT{u')du\ (36a) 



gtiu) = '-^m- (36b) 

Here, 6{x) is the standard step function, taking the value 1 or according to whether x is 
positive or negative, respectively. Equations (^) and ( p6D describe how the function ip^^Q 
can be constructed given any form r(-u) for the compact initial outgoing pulse. 

Eqution (^) implies that the wave V'i=o sharply "cut off" at retarded time u = 0. 
This effect is due to ingoing and outgoint waves destructively interfering with each other 
at the origin of coordinates (r = 0) inside the fiat region internal to the material layer, at 
retarded times m > 0. [For a scalar field in the "complete" Schwarzschild model we found 
in Ref. [§ an exponential decay of the waves at m > rather than a sharp cutoff — thus 
the main support of ipn^Q remains compact even in this more sophisticated (and much more 
complicated) model.] 

Obviously, since V^i=o strictly compact in retarded time (it is supported only in the 
range mq < m < 0), it does not contribute to the overall late time (large u) radiation at 
null infinity. Rather, the function ip^^f^ will serve [via Eq. (|5D|)] as a source to higher-order 
(n > 1) terms of the iteration scheme, which will form the late time tail of decay, as we 
show below. 



C. Green's function of the iteration scheme 

Using the Green's function method, we formally have, for each of the functions ipn^^Q, 

iIjI{v,u)= r dv' r du'G\v,u]v',u')Sl{v',u'), (37) 

Jo J uq 

which allows one, in principle, to calculate these functions one by one, in an inductive 
manner. Here, G{v,u;v',u') is the time domain Green's function in Minkowski spacetime, 
defined as satisfying the equation 

G[uv + V^ir)G' + Ro{r)G\, = 5{u - u')5{v - v'), (38) 
with the causality condition 

G\u < u') = G'{v < v') = 0. (39) 

To solve for G, we use a straightforward generalization of the method used in Ref. [0. 
First one shows that for Eq. (|38| ) to be consistent with the causality conditions (|39|), one 
must have G{v = v') = [{v' — u)/{v' — u')Y and G{u = u') = [{v' — u')/{v — u')^. This 
establishes a characteristic initial-value problem for the Green's function ai u > u' and 
V > v'. Then, with the help of Eq. (|36|), one can obtain (see [0 for more details) 

G^v, u- v', u') = E Af ^. ^"'^;7 e{u - u') e{v - v') eW - «), (40) 
j=o [y -u) 
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where the function g{u; v', u') is given by 



g{u-v\u') = ^ , ,M ' ^1 

(f + s)! [v' — u'Y 

and where the j differentiations of this function are with respect to u. The factor 6{v' — u) in 
Eq. (^) is related, again, to the presence of the origin of coordinates inside the flat internal 
region. 

Finally, a general comment should be made about the application of the iterative scheme 
to the Kerr spacetime: Here (unlike in the Schwarzschild case) the characteristic surfaces of 
the iterative scheme (i.e. the surfaces of constant u or constant v, which are not null but 
rather timelike) do not coincide with the actual characteristics of the Kerr geometry. Thus, 
strictly speaking, the "causality condition" stated in Eq. (|39|) does not hold for the actual 
Green's function in the Kerr spacetime. However, the surfaces u =const and v =const do 
approach the actual null characteristics of the Kerr background at large distances. Recalling 
that the form of the late time radiation at null infinity is shaped mainly during the propa- 
gation at large distances, it is reasonable to expect that the above problem will not crucially 
affect the validity of our results. 



D. First-order iteration term 

The first contribution to the late time (large-w) tail at null infinity comes from the 
function ipn^i- This contribution also turns out to be the most dominant one, with those 
of the functions ip^yi smaller by one or more factors of M/\uo\ (recall that in our model we 
have |mo| ^ M). It is therefore of special importance to analyze in detail the behavior of 
ipn=i, and derive its late time form at null infinity, as we shall do now. 

The function V'i=i is calculated from Eq. (|37|), with n = 1. Since we only look for the 
behavior at null infinity, we take the limit f ^ oo of this equation. At that limit, the Green's 
function appearing in the integrand is dominated by merely the j = I — s term of the sum 
in Eq. (1^) : 

G\v ^ oo) = [g{u- u\ v')f-'^ e{u - u') e{v' - u) (42) 

[I — s). 

(where the derivatives of g are with respect to u). 

With Eq. (^) we now obtain for il'L=i ^^^^ infinity (for each /), 



{v' — uy^'^{u — u 



(Is) 



'4^n=iiu,v ^ oo) = du dv'^ — ^ Sl^^{u\v') (43) 

(where, again, the / — s derivatives are with respect to u). In this expression, the lower 
limit of the integration over v' was set io v' = u due to the factor 6{v' — u) appearing in the 
Green's function, and the upper limit of the u' integration was set to m' = in view of the 
compactness of ipn=Q (implying the compactness of the source 5*^^^ as well). 

The source function 5*^^^ is calculated from Eq. (pID, with n = 1. It contains, in general, 
contributions from the modes = /, / ± 1, / ± 2 of ipn=Q- These contributions to V'i=i [via Eq. 
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(1^)] are additive, and may be calculated one by one. The details of this calculation are given 
in the Appendix. In brief, it contains two steps for each of the above contributions: First, 
the definite integration over v' is carried out explicitly. The integrand of the remaining 
u' integration then becomes a finite sum over derivatives of g^lu'), each multiplied by a 
power of {u — u'). In the second and final step we use successive integrations by parts to 
eliminate the derivatives of {u'), with all resulting surface terms vanishing in virtue of the 
compactness of this function. (This procedure is clarified in the Appendix.) 

The following is a description of the outcome from the above calculation. Let us denote 
by ip^Z^j the contribution to the mode / of ipn=i at null infinity from all terms in Sl^^^ 
associated with the mode I' of ipn=o- We then find at u ^ — ^O) to leading order in M/u 
and in uq/u (see the Appendix for details), 

where the a's and /3's are constant coefficients (depending on s,l,m), and where /q is, for 
each I, a simple functional of Qq: 




For a scalar field, s = 0, the contributions 4'n=i^ vanish (as the coefficients a^.]^ are propor- 
tional to s). Also, there is no contribution ipl^^'' for Z < /q + 1 and no contribution 'ipl^-[^'' 
for / < Zo + 2, where Iq is the lowest radiatable multipole mode for given s and m: 

Iq = max(|s|, \m\). (46) 

Equation (^4]) implies that for any given s and m, the most dominant mode of ipn=i at 
null infinity, at large u, is Iq. It also tells us that the dominant contribution to this mode 
comes solely from this mode itself (namely, it is tp^Zl''^), since in this case there are no 
contributions from lower modes, and the ones from higher modes are negligible at large u. 
Hence, the decay of the dominant mode of ipn=i at null infinity (for any given s and m) is 
described at large retarded time u by 

^n=i{u, v-^oo) = (47) 

to leading order in M/u and in uq/u. Recall that the function g!o{u) and, hence, the 
functional Jq" are nonvanishingj^ as long as we make the assumption that the initial data 
have a generic angular form (so that it particularly contains the lowest multipole mode Iq for 
any given s and /). The case where the lowest modes are initially missing will be discussed 
in brief below. 



1 + Pli 1^ (^) 



(44) 



° This is true unless the initial data are very finely tuned such as to make the integral in Eq. (45) 
vanish. 
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E. Late time tail at null infinity 

Equation ( pTf ) provides a formal means for calculating the higher iteration terms, ipn>2- 
However, exact analytic calculations become very tedious already for the n = 2 term. In 
the case of a scalar field in Schwarzschild spacetime we have explicitly derived V'n=2 
null infinity, and showed that it exhibits the same power-law decay as i'n^i at large u, yet 
with an amplitude reduced by a factor proportional to (M/|mo|) <C 1. For this case, analytic 
considerations suggested that a similar reduction of the amplitude by a factor oc M/uq 
occurs also for n > 2, whenever n is increased by 1. This conclusion was verified numerically 
for the first few iterative terms [0. Our numerical calculations also indicated that the sum 
of iterative terms ipn seems to converge rather fast at null infinity for large \uq\/M (say, in 
the order of 100). 

We now proceed under the assumption that the same considerations also apply in our 
case, of general-s fields on the Kerr background. That is, we assume that (for any given 
s,/,m) the functions 'ipn>2 decay at null infinity at large u with the same tail as V'n=ij Y^t 
the amplitude of these functions is smaller by at least one factor of order 0{M/uq). This 
assumption seems plausible, because the above property of the iterative scheme [namely, 
the scaling of tpn as oc (M/uq)"] seems to stem from the basic structure of the iteration 
procedure, rather than from the details of the source function Sn, which distinguishes the 
Schwarzschild, s = 0, case from the more complicated case studied in the current paper. 

Adopting the above assumption, we conclude that, for large \uo\/M, the "overall" func- 
tion -0' is well approximated at null infinity at large u by merely the term ip^^i- In particular, 
lo is the dominant mode of ip there, for any given s and m. By virtue of Eqs. (P) and (^T]) 
we then finally obtain for the Newman- Penrose field \E'''™ (for any given s,m), 

^sm ^ ^lo jio ys/om(^^^ ^) ^-2s-i ^-{io-s+2) (^jominaut mode at scri+, large u). (48) 

This result is accurate to leading order in M/u, in Uq/u, and in M/uq. In the generic case 
where the initial pulse includes all values of m, we find that the behavior is dominated by 
the modes with < \m\ < \s\ and / = which decay at null infinity with the late time tail 
for s > or m^"*"^ for s < 0. 

One may also ask about the behavior of the other, faster decaying, modes at null infinity. 
From Eq. (|4^) we find that (for any given s and m) the modes I > Iq of ipi are also "fed" 
by strong contributions coming from modes of smaller I: In general, the function has 
leading-order contributions from the modes / — 1 and / — 2 of il'n=Q- In the exceptional scalar 
field case (s = 0) this contribution comes only from the / — 2 mode (provided / > /q + 2), 
and provides the same tail as the contribution from the mode / itseld, namely u^^^^^^'^\ We 
thus have, for all / > /q in the scalar field case, V^^^i oc ^"('"■^+2) at null infinity, large u. 
Under the above assumption that the "overall" field is well approximated there by ipi 
(for all /), we conclude that in the scalar field case, the decay of any of the modes I >Iq aX 
null infinity is given at large u by 

where F'™ are the spherical harmonics. 

As to the non-dominating modes of the s 7^ fields: Equation (0) suggests that 
these modes would exhibit not a strict power-law tail but rather a tail of the form 
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(X u {!■ s+2) X ln( ^ )■ We feel, however, that this result cannot be taken as conclu- 
sive, and needs a further support (e.g. from numerical analysis). (We comment that such 
logarithmic dependence does not arise from the frequency-domain analysis in Ref. ||2^.) We 
emphasize our conclusion, Eq. (|48|), that the leading-order tail at null infinity, belonging to 
the most dominant mode, decays with a strict power-law. 



IV. LATE TIME EXPANSION 



The target of this work is to explore the behavior of the fields at late time anywhere 
outside the KBH (and along its EH). To that end we shall apply the late time expansion 
scheme, a version of which was used in Ref. |§] to analyze a scalar field in the Schwarzschild 
case. 

We assume that at late time, the fields admit an expansion of the form 



vl>^(t;,r,^,(^) = 5: 



fc=0 



l=\s\ m=—l 



slm I 



Im I 



-ko—k 



(50) 



to which we shall refer as the Late Time Expansion (LTE). Here, ko is a constant parameter 
which we later determine.^ As we show in this paper, the LTE is consistent with the field 
equations, with the tail form at null infinity, and with regularity requirements at the EH. 
We adopt an expansion in 1/v, rather than in 1/t, because it appears to be more adequate 
for analyzing the behavior near and along the EH (as the coordinate v, unlike t, is regular 
through the EH). 

Inserting the form (|50D into the master field equation (|^), and collecting terms of common 
V powers and of common multipole numbers / [with the aid of Eqs. (|12|) and (p!^], yields an 
ordinary equation for each of the unknown functions 



slm I 



D 



slm 



-^slm I 



7slm 



in which Z}*'™ is a differential operator given by 

a^im? + 2isma{r — M) 



D' = A^ + 2{s + l){r-M)^ + 
dr^ dr 



and the source term Z\. reads 

Zi = 2{k, + k-l)[{r^ + a')'^ + 



dr 



A 



2M[s(r2 



-{l-s){l + s + l] 



(51) 



(52) 



imar\ 



A 



lasc, 



F 



k-1 



las 

.2 



c^-Fi-\ 



+ y(A:o 



{c'^Ftl 

k - 2) {cl^Ft 



k-2 



+ Ci^Fitl + CiFitl + CiFtl + 



" + k-2 



k-2 



1-2 
k-2 



(53) 



' It should be noted that hy definition the parameter /cq does not depend on /,m: whereas in the 
Schwarzschild case a separate parameter has been defined for each mode /,m, in the present 
paper a single parameter /cq is related with the overall field 
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with -F^<o = 0) with the various coefficients and given in Eqs. (|T3D and ([T5|). 

An essential feature of Eq. ( ^l]) is the fact that it actually constitutes a hierarchy of 
effectively decoupled equations, as each of the functions Fk^Q satisfies an inhomogeneous 
equation whose source depends only on the functions -Ffc'<fc preceding it in the hierarchy. 
The ffist function, F^^q, obeys a closed homogeneous equation. Thus, in principle, we may 
solve for all modes of all functions Fk one by one, starting with k = 0. For each k, one 
should be able to solve for all modes I of Fk, and then carry on to + 1. 

Now, Eq. (|5T|) is a second-order differential equation for each of the various functions F^. 
In principle, to determine these functions, proper boundary conditions should be specified at 
the EH (which is, mathematically speaking, a "regular singular point" of the equation) and 
at space-like infinity. (These two boundary conditions should determine the two arbitrary 
parameters which occur in the general solution for each of the function F^.) The behavior 
of \E''' at infinity is known from the previous section, and in Sec. [VII| below we discuss the 
matching to this asymptotic region. At the EH, the only obvious requirement concerns 
the regularity of the physical fields there. In the rest of the present section we obtain the 
boundary conditions for the functions Fj;. at the EH, based on local regularity considerations. 



Boundary conditions at the EH 

One expects "measurable" physical quantities to maintain a perfectly regular behavior at 
the EH, which is a surface of a perfectly regular local geometry. Accordingly, the components 
of the Weyl and Maxwell tensors should be perfectly smooth through the EH, provided these 
components are expressed in a coordinate system regular at the EH. To construct boundary 
conditions for the scalars at the EH, it then remains to relate these scalars to the regular 
components of the Weyl and Maxwell tensors. 

To that end we must ffist write the tetrad basis (^), used to construct the scalars \E''*, in 
EH-regular coordinates. Recalling that the BL coordinates t and ip go irregular at the EH, 
we introduce the Kruskal-like coordinates 

U = -e-"+" and V = e^+\ (54) 

and the regularized azimuthal coordinate 

ip+ = (f — Q+t (55) 

(see Sec. 58 in Ref. [0]), where 



+ 2Mr+ ^ ^ 

(f2_i_ is the "angular rate of inertial frame dragging" at the EH.) In the EH-regular coordinate 
system (V, U, 9, (p+), the components of the ingoing and outgoing tetrad legs have the EH- 
asymptotic forms 

I, A-ie''+-[l,0,0,0], . . 

n'' oc A e-'"+^ [0,1,0,0]. ^ ^ 

Recall that v is regular at the EH, and that A = there. 
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Now, the construction of the scalars involves |s| projections of the Weyl and Maxwell 
tensors on the tetrad legs l'^ (for s > 0) or n'^ (for s < 0). Since the components of these 
tensors must be EH- regular in the coordinate system {V,U,6,ip^), then by virtue of Eq. 
(^) we find that A'*\l/* must be EH-regular as well. 

To formulate the regularity condition for the functions Fl{r), we note that the functions 
Y'^^'^iO, ip) in Eq. (|50D are irregular at the EH for m 7^ 0, due to the factor e*™'^: We have 



g-imn+r.^ (58) 



where the factor in square brackets is EH-regular, but the following factor oscillates rapidly 
towards the EH. From Eq. (^) we thus find that it is the quantity /\«e-«'"f^+^*^^ which 
must be perfectly regular at the EH (for all k). 

Hence, the regularity condition at the EH can be phrased as follows: Define the "physical" 
variables 

^' = A"^' and = A"e-*'"^+"*F^. (59) 

Then, we must have that (for all k) 

'^'^ and are smooth functions at the EH. (60) 

Mathematically, we shall require these functions and all their derivatives with respect to an 
EH-regular coordinate (such as r or U) to be continuous through the EH. 

Equation (pO| ) constitutes the required boundary condition at the EH for all functions 
Fl{r). 



V. GLOBAL SOLUTION FOR A; = 

The dominant late time decay at worldlines of fixed r is described by the = term 
of the LTE, Eq. (0).0 In this section we derive an exact analytic expression for this term 
(namely, for all modes / of the function F^^q). Since Fj^^Q is a solution of a second-order 
differential equation, it shall contain two arbitrary parameters. One of these parameters will 
be determined in this section by the regularity condition at the EH. The other parameter 
will be determined in Sec. [VI 1| through matching at infinity. 



By definition, Zl^^ = 0, and the function F^^q admits the homogeneous equation 

FLo = 0, (61) 



with the operator D given in Eq. (|5^) . This is nothing but the static field equation in Kerr 
spacetime. The static solutions play an important role in our analysis, for two reasons: (i) 
As just mentioned, the late time behavior can be approximated by knowing -F^^O' which 



^ Actually, in Sec. VI we discuss an exception to this statement: For m = modes of s > fields, 
the behavior along the EH is dominated by the k = 1 term. However, in this case too, the decay 
along Unes of constant r > r+ is still dominated by the A; = term. 
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must be a static solution, and (ii) we shall use a basis of static solutions in constructing the 
functions FI^q using the Wronskian method. 

For reasons that will become clear below, we continue by treating separately the cases 
s 7^ and s = 0. Also, for s 7^ we will consider separately the cases m 7^ (nonaxially 
symmetric modes) and m = (axially symmetric modes). 

A. s 7^ fields: Nonaxially symmetric modes 

For s 7^ and m 7^ 0, a basis of exact solutions to the homogeneous static field equation 
(|6T|) is given by 

(f)r = {z+/z_y"''^ F{-1 - s,l - s + l;l - s + 2im-f; -z+), (62a) 

(j)ir = {z+Zz.)''""^ F{-1 + s, / + s + 1; 1 + s - 2im7; -z+), (62b) 
where z± are the dimensionless radial variables defined in Eq. ([21|), 



7 = = 7, — . (63) 

r+ — r_ 

and F denotes the hyper geometric function (We use this notation because, as we discuss 
below, (pr is physically regular at the EH, whereas ipir is irregular there.) 

The hypergeometric function F(a, 6; c; y) (where a, 6, and c are complex parameters and 
7/ is a complex independent variable) admits the series expansion 

F(a,6;c;y) = l + f: (64) 



n=l n\ 



(see, e.g.. Sec. 2.1.1 in ||5^), where 

(a)„ = a(d + 1) ■ ■ ■ (a + 72 — 1) = V{a + ri)/V{a) (65) 

is the "rising factorial" .0 Two results arising from Eq. (|6^ are that (i) the hypergeometric 
function is not defined if c is a non-positive integer [as in this case a zero factor occurs in 
the denominator in Eq. ([6^)], and that (ii) if either d or 6 are non-positive integers, the 
expansion ( |6^ terminates, and the hypergeometric function becomes a polynomial of order 
—a or —6, respectively. 

Item (i) above implies that for m 7^ both solutions ( |62a[ ), ( |62b| ) are defined; however, in 
the case m = (which we treat separately below) only one of these solutions is defined (^^ 
for s < 0, or (pi^ for s > 0). We further find, by item (ii) above, that both hypergeometric 
functions in the solutions ( |64D are simply polynomials of z+ (and thus of r too). 



^ Note that {a)n is well defined even when the expression involving the gamma functions is not. 
In this case, of course, the second equality in Eq. (|65|) is invalid. 
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For m 7^ 0, the general static solution is constructed from the two basis functions (|62D. 
With the help of Eqs. (^21) , (|6^), and (^), one finds the asymptotic forms of these functions 
to be 

s^imn+r, ^ _^ _ 0), ^gg^ 




as r, z I ^ OO, 



and 




as ^ -OO (A, 2+ -> 0), ^^^^ 



r' as r, 2;+ ^ 00 



where the coefficient B^"^ is given by 



,i„ (2i)!r(l-s + 2,m7) 



We now use Eq. (|9]) to construct the "physical" fields 4>r = A'^e ^"^^+^*(j)r and (^jr = 
^Sg-jmO+r,0,^ associated with the static solutions (pr and (pir. From the EH-regularity 



criterion (|60D we then learn that the static solution 0^ is physically regular at the EH, 
whereas (pir is irregular there: For s < the irregularity of (pir is obvious, as it diverges like 
~ A^l'^l at the EH. For s > 0, (pir is continuous through the EH, yet its sth derivative with 
respect to U (which is a regular coordinate through the EH) diverges there like oc e~^*''"^+''* . 
Higher order U derivatives of (pir are unbounded in magnitude at the EH. 

Obviously, (pr is the only solution of the homogeneous static equation (up to some global 
constant factor) which is physically regular at the EH, because any combination a(pr + b(pir 
with 6 7^ will be irregular there. In a physical setup where a static source presents outside 
the BH, the field must behave as oc (pr near and through the EH. In vacuum, (pr is the only 
global static solution physically regular at the EH. This field does not vanish at infinity 
(where it behaves like oc r'"''); hence it cannot represent a physical static perturbation. 
There exists a static solution which dies off fast enough at infinity [this is the solution 
{B'-_!!l)*(pr — B''J^(pir, which dies off as oc r^'^^^^ at large r [37|], yet this solution is physically 



irregular at the EH. Similar results apply also to the static m = modes, which we study 
below. We thus conclude that there cannot exist physical vacuum static modes outside the 
Kerr BH. This, of course, is a manifestation of the "no hair" principle. 

In the framework of the LTE, each of the functions Fl{r) must be subject to the EH- 
regularity criterion (|5UD. Since F^^q must be a static solution, it thus have to be proportional 
to (pi. To conclude the above discussion we therefore take 



r' 



Fi^,{r) = a[(t>[{r). (69) 
The constant Oq is to to be determined in Sec. |VII| by matching at null infinity. 



B. s 7^ fields: Axially symmetric modes 

We remarked earlier that in the case m = only one of the basis functions ( |62aD , 
( |62b| ) is defined — the one in which, depending on the sign of s, the third parameter of the 
hypergeometric function is a positive integer. Denoting this function by 0^^°, we have 
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^rn=o ^ ^s,l + s + l;s + l] = </)+ for s > (70a) 

and 

0™=° = A-'F{-1 -s,l-s + l;-s + l; -z+) = 0^ for s < 0. (70b) 
The asymptotic forms of this solution are 

C-M ^ { f-. : > »j at the EH (71) 



and, for both s > and s < 0, 



m=Oi'^\ ~ r)l,m=0 l~s 



[r 



as r, ^4. — > oo. (72) 



Recall that F is simply a polynomial of 2+ (and of r), and thus so is 0™^°. From Eq. (ffOD 
we find that the "physical" field (/)™^° = A'*^^^'' is also a polynomial, and therefore, clearly, 
0™^° is physically regular at the EH. 

We still have to construct a second independent basis static solution for the m = case. 
(This will allow us to tell whether or not the above solution 0™^° is the only regular one.) 
Fortunately, at that point we can benefit from the work already done in Ref. for the 
Schwarzschild case: When expressed in terms of the variable z+ (rather than r), the static 
field equation (|6TD for m = takes exactly the same form as for the Schwarzschild BH [see Eq. 
(21) in [^], where in the latter case we use the variable z = (r — 2M) / (2M). Therefore, each 
static solution in the Schwarzschild spacetime becomes a static axially symmetric (m = 0) 
solution in Kerr spactime, upon replacing z —>■ z^. Moreover, in terms of the variables z (in 
the SBH case) and z^ (in the KBH case), the EH- regularity criterion becomes the same for 
both spacetimes, and thus the classification of regular and irregular solutions at the EH is 
also conserved. 



As a second basis function we then take the static solution given in Eq. (24) of Ref. p8|: 



im=0 
ir 



~Ais zl' zZ^-^Fil - s + 1, / + 1; 2/ + 2; zZ^] 



for s > 



(73a) 



and 



I m,=0 
ir 



Ais zZ^-'-^Fil + s + 1, / + 1; 2/ + 2; zZ^) = (p'r for s < 0, 



(73b) 



where we have replaced z with Zj^ (and thus z + 1 with z_). Here, Agi is a normalization 
factor, 



Asi = 1/F(/- |s| + 1,/ + 1;2/ + 2;1) 



/!(/ + |s|)! 



(74) 



(2/ + l)!(|s|-l)! 

[cf. Eq. (46) in Sec. 2.8 of ^8[], chosen such that 0™""° takes a simple asymptotic form at 
the EH (see below). From Ref. [^] we also know that 0^^° admits the following series 
expansion near the EH: 
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A-^(l + a+A + ---)+/3 0+1112+ (fors>0), , . 

{l + a-A + ---)+(3(f);\nz+ (for s < 0), ^ ^ 



in which the coefficient f3 is nonvanishing, 

3. (zir™^ pe) 



It is clear (e.g. by comparing the EH-asymptotic forms) that the two solutions 0^^° and 
0™""° are independent, and thus form a complete basis of solutions. As already explained in 
Ref. p8|], 0[J^"° is physically irregular at the EH: For s < the "physical" field 0™""° diverges 



there as A~I*L For s > 0^^° is continuous through the EH, yet its sth derivative with 
respect to U diverges there (as oc lnz+). Therefore, Eq. ( |69D applies to m = as well, where 
for this case the function 0^ is given by Eq. (|7^). 

It is instructive to compare between the asymptotic behavior of the s 7^ static solutions 
in the case m 7^ [Eqs. (|66D and (|67D ] and in the case m = [Eqs. (^) and (^)]. Focusing 
on the asymptotic dependence on A (and ignoring for this discussion the oscillatory factor 
presents in the m 7^ case), we find that for the s > fields the regular and irregular 
solutions "switch roles" : For m 7^ modes the regular solution is the one that behaves like 
A~'^ at the EH and the irregular solution is the one that behaves like oc const there, whereas 
for m = modes the opposite is true. Such an interchange of roles does not occur in the 
case s < 0. This effect is explored and explained in detail in Ref. 



C. Scalar field case (s = 0) 

For s = we use the new radial variable 

2r - r+ - r_ 

z = z+ + = (77) 

r+ — r_ 

(note the relation z = 2z+ + 1 = 2z- — 1), to write the static field equation (|6T|) in the form 

(1 - z^)F"{z) - 2zF'{z) + [4m^^\l - z^)'^ + /(/ + 1)] F{z) = 0, (78) 

where a prime denotes d/dz. This is the familiar Legendre's differential equation (see, for 
example. Sec. 3.2 in Ref. |3^). Two independent solutions to this equation are 

0^=0 = {z+/z_Y'^^ F{-1, l + + 2imr, -z+) (79a) 

and 

is=0 T?/l /o I A. I 1 1 /O I A. I 1 /O. 7 I Q /O. ;;;-2\ 



— F{l/2 + im-f + 1, 1/2 + + 1/2; I + 3/2; z'^), (79b) 



which are (up to a customary normalization) the associated Legendre functions of the first 
and second kinds, P^iz) and Qf (^), respectively, with /x = 2im7.[^ 



^'^ Note that the two independent solutions 0^^^^ and 0^^°, Eqs. ( |62aD and (|62b| ), degenerate to 
the single solution ( [79 al) in the scalar field case (i.e., when setting s = 0). 
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Using the EH-asymptotic relations z+ ~ gSs+r, _ gf2+r,/7 ^^^^^ ^_ ~ 1, and recalling that 
the hypergeometric function appearing in Eq. ( [TQap is a polynomial (of order /) of we 
find - 6*™^+''* near the EH. The "physical" field 0^=° = e-*™^+^*0^=° associated with 
this solution is therefore regular at the EH. 

It remains to verify that 0*^° represents the only physically regular solution for the scalar 
field (up to a constant factor). For m 7^ 0, the solution (pl""^ and its complex conjugate 
(0^^")* constitute a complete basis of static solutions. The "physical" field e~*'^^+''*(0^^°)* 
becomes indefinite at the EH (where it behaves like g-2«"if^+»'*^^ g^^^j ^f^Q that for m 7^ 0, 
0^^° is indeed the only physically regular solution. In the case m = 0, 0^^° becomes 
real (it is then the Legendre polynomial, up to a normalization), and a complete basis of 
solutions is given by (0*"°,0|^°). These two basis functions then admit a relation of the 
form 0^jr°(z) oc 0*"°(z) X ln(2+/2;_) + polynomial in z [see Eq. (24) in Sec. 3.6.2 of [^] . 
Since 0^"° is physically regular, it is therefore clear that is physically irregular. 

In conclusion, because the function F^'Zq must be a static solution physically regular at 



the EH, it must be proportional to 0* °. Therefore, Eq. ( p5| ) applies to the scalar field too, 
with 0r given in Eq. (|79a|) . 



VI. LATE TIME BEHAVIOR AT THE EH 

The LTE, Eq. (|50|), is an expansion in inverse powers of advanced time f, with r- 
dependent coefficients. Since along the EH itself r is constant and v takes finite values, 
this expansion seems especially convenient for analyzing the "late time", v ^ M, behavior 
of the fields at the EH. Potentially, this behavior should be described by the A; = of the 
LTE. However, a possible divergence or vanishing of various "coefficient" functions Fk{r) 
at r = r+ may alter this simple picture, leading to a different prediction for the late time 
power-law decay at the EH. Indeed, as it turns out in this section, there is a case (the one 
of s > 0, m = 0) in which the term = 1 is found to dominate the term = at the EH. 

It is therefore important to analyze also the behavior of the k > 1 terms at the EH. 
This task is further motivated by our wish to verify that the LTE is fully consistent with 
regularity at the EH: It will be shown that for each k there exists a solution Fk physically 
regular at the EH. These EH-regular functions will then construct, via the LTE, a field 
representing a physical perturbation which is regular along the EH at all v. 

With the above motivations in mind, we first derive in this section expressions for the 
EH-asymptotic behavior of each of the functions Fk>i. The results are summarized in Table 
I below. We then discuss the application of these results to the late time tail along the EH. 



A. Behavior of the A; > terms at the EH 

Each of the functions F^-^q admits the inhomogeneous equation D^Fl = Zj. [Eq. (pl|)], 
and is subject to the EH-regularity condition (|60D. For each A; > 0, the general solution to 
Eq. ( pID has the form 

Fl>oir) = dkMr) + bkMr) + <pf{r), (80) 
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where and bk are (yet) arbitrary coefficients, (f>r and (pir are two independent homo- 
geneous solutions (those derived in the previous section), and 4>]l^{r) is a solution to the 
inhomogeneous equation. 

For each /c > 0, an inhomogeneous solution cp]^ is given by 

0, (r) = 0,(r) j dr - 0,,(r) j dr , (81) 



in which 



W = A-'-^ (82) 



is the Wronskian associated with the homogeneous equation D'-Fl = 0. We can now make 
use of the relation 



cpir = -Mr) J ' (p;^{r')W{r')dr' (83) 
to re-express the above inhomogeneous solution in a more convenient form, as 



lih I 



dr' rdr" ^^^'^'^^■^'"^ ^ ^ (M) 
'r, [(pr{r')Y W{r") A{r") ' ^ > 



where ri and r2 are constant integration limits.|^ This form is obtained from Eq. (|8T|) by 
first substituting for (pi^, using Eq. (^), and then integrating the resulting expression by 
parts. It is advantageous in that it only involves the homogeneous solution 0^, which is of 
a more simple form than in all the cases considered in the previous section. 

Equation ( p^ can be used, in principle, to calculate all functions Fl{r) in an inductive 
manner. In general, for each > the source function depends on various / modes of 
the functions Fk-i and Fk-2i which are to be calculated at previous steps of the induction 
procedure. As we show below, for each /c, the value of one of the coefficients or is 
dictated by regularity at the EH. The other coefficient is to be specified by matching at 



large distance, as we explain in Sec. VII. 



We now use Eq. (^) to obtain the EH-asymptotic forms of all functions F^^q. The 
special case s > 0, m = shall be treated separately from all other cases. 



1. The case s < and the case s > 0, m ^ 

We start with k = 1. The source ^^=1 is calculated from Eq. (^), in which, according 
to the results of the previous section, we set F^^q = a^cpl. For s < 0, and also for s > 
with m 7^ 0, we find by Eqs. ( |62a|) , ( [70 b|) , and (|79a|) that the function FI^q has the form 
-F^^Q oc 0', = A"^"^**^"^ X fo{r), where /o('") is a certain function analytic at the EH and 
nonvanishing there. Substituting this form in Eq. ( ^3]) yields 



The constant limits ri and r2 will be specified as convenient for each of the various cases 
analyzed below in separate. Of course, changing these limits amounts to adding a homogeneous 
solution to Fj^, which is merely equivalent to re-defining the coefficients or bk in Eq. (|80| ) 
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ZU{r) = A-+^-^ X Mr), (85) 



where fi{r) is a function analytic at the EH.0 

With Eqs. (|8|) and (H), Eq. (H becomes (for k = 1) 



^t. = Mr) I *' I ||(;,)|.]|^(,,)|.,. X /.(r"). (86) 

where fi{r") is analytic at the EH, and where we have specified the lower limit of the 
integration over r' as ri = r+. Integrating over r" [recalling A = {r ~rj^.){r ~r_)], we obtain 



r 



£ rfr' (g + c [A{r')Y-'-'^^^) X 7i(r'), (87) 



where /i(r') and fi{r') are analytic at the EH, c is a certain nonvanishing constant, and c 
is an integration constant, which, of course, depends on the value of the lower integration 
limit r2. In the case s < 0, a convenient choice is r2 = r+, which makes c vanish. For s > 
we take r2 = 2r_|_ (say), as the choice r2 = r+ is forbidden^. In that case, the contribution 
proportional to c to the integral in Eq. ( |57D must coincide (up to a multiplicative constant) 
with one of the static solutions, (pr or because changing the integration limit r2 (thus 
changing c) amounts to adding a static solution to Eq. (|80|) . By integrating over r' we find 
that this contribution has the asymptotic form oc cA"*™''' oc ce~*"*'^+'^* at the EH; hence 
it must admit the global form oc c(j)ir [see Eq. (0)]. The term proportional to c in (/'^'Li 
can therefore be absorbed in the term fei^jr of Eq. (0), by re-defining the coefficient hi. 
One is left with the contribution proportional to c, which, after integrating over r', reads 
0r ■ A fi{r), where fi{r) is a function analytic at the EH. 

From Eq. (|80D we now find that Fl^^ = ai(j)r + &i0ir + 4>r^fi{f)- The EH-regularity 
criterion (|60|) then dictates bi = 0, finally leading to 

FLi=0.[«i + A-/i(r)]. (88) 

Recall (i) that /i(r) is analytic at the EH, and (ii) that the regular solution is physically 
regular there in the sense discussed in previous sections. This implies that the function F^.^^ 
of Eq. (^) is physically regular at the EH (namely, the "physical" function Fl^i associated 
with is mathematically regular there). Note also that, as far as the leading order term 



In deriving Eq. ( |85| ) one should notice that when F^^q is substituted in Eq. (|53|), the term 
containing dFl^^/dr and the one containing F^^q/A cancel out at the leading order in A. As 
a consequence, leading-order contributions to Z^.^^ arise from all terms in ([5^), including the 
interaction terms. 

In the case s = +2 we can make c vanish by taking r2 = oo. However, for s = +1 no choice of 
r2 nullifies c. 
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in A is concerned, Fl^^ has the same asymptotic behavior as FI^q at the EH, which is that 
of the regular static solution: oc A~'^e*'"^+'"*. 

We can now carry on in an inductive manner, and analyze the terms k > 2. In general, 
let us assume that for a given k' > 2 we have for all k < k' 

Fi = (j)r[ak + A-Mr)], (89) 

where fki''^) are functions analytic at the EH. (This form was already verified above for k = 
and k = 1.) By substituting in Eq. (|53D it is straightforward to show that 

Zi,(r) = A-^+^-^x/fc,(r), (90) 

where /fc'(r) is analytic at the EH. This, following the same calculation as for leads 
to Fl, = ak'ipr + bk'(pir + 4'rAfk'ij'). The EH-regularity condition (160|) then dictates hk' = 0, 
and one finds that Eq. ( p9D is also valid for k = k'. Thus, by induction, Eq. ( pQ]) is verified 
for all k>0. 

In conclusion, we have constructed EH- regular solutions for each of the functions Fl{r). 
It was found that for all modes /, all functions Fl{r) behave near the EH like a regular static 
solution, oc ipr', namely, they all admit the asymptotic form 

Fl{r) = ak A-'e*'"^+"* near the EH, for all A; > 0. (91) 



2. The case s > 0, m = 

In this case, the function FI^q = ao0^ is simply a polynomial, admitting the EH- 
asymptotic form = aoA" [see Eq. (|7lD]. To obtain the source Zl^^, insert -F^=o 

into Eq. (^). This yields 

^s>o,i,m=o ^ ^Mskoaoirl - a^)A~' + h{r), (92) 

where h{r) is a function analytic at the EH. Note that now, since -^^=0 ^ polynomial, the 
asymptotic form of S]^^^ at the EH is dominated by merely the term oc A~^F^^q in Eq. (^3|), 
while the other terms (including the derivative term and the interaction terms) contribute 
only to higher orders in A. 

To calculate 0^=1; "we substitute S\^^^ in Eq. (P^. Recalling that in the present case the 
function (pr is a polynomial, one finds 

= AMsk,a,{rl - a^mr) [[ dr' £ dr''^-^^^^^^^ (93) 

where h{r") is analytic at the EH. (Here, we have specified one of the integration limits, 
r2 = r+.) Carrying the integration we arrive at 

2 2 

= 4Mkoao ~" 0+(r) T dr' \A-\r') + hi{r') 
r+ — r_ Jri L 

= ji(f)+ In z+ + hi{r), (94) 
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in which hi{r') and hi{r) are analytic at the EH, and where the coefficient 71 is given by 



71 = AMkoao + 0. (95) 

(r_|_ — r_Y 

By virtue of Eq. (^) we finally obtain for F\.^^ (in the case s > 0, m = 0), 

Fi=i(r) = a+0+ + 6+0+ + 0ti 

= a+0+ + hX\^-\\ + a+A + ■■■)+ In z+] + 71 0+ In 2+ + /ii(r). (96) 

Now, the EH-regularity criterion (|60D forces the "physical" function F\^^ = A^Fj^^^ to 
be perfectly smooth at the EH (where = 0). This implies that must contain no 
logarithmic terms of the form oc lnz+: If such a logarithmic term is present, Fj^^i would 
indeed be continuous at the EH, yet its s'th derivative with respect to r or U (which are 
regular coordinates at the EH) would diverge there. We therefore find that the regularity 
condition dictates the value of the coefficient bf: 

bt = -11/ ~P + 0. (97) 
Hence, from Eq. (^) we obtain the form 

FLi(r) = A-%(r), (98) 
where h\[r^ is a function analytic at the EH, satisfying 

/ii(r = r+) = 6+ ^ 0. (99) 



[The analytic function /ii(r) contains also the terms 01"^^ and /ii(r) appearing in Eq. 
multiplied by A^. Note that the polynomial homogeneous solution does not affect the leading 
order term of F\,^^ at the EH.] It is important that /ii(r) does not vanish at r = r+: It implies 
that at the EH itself (in the case of s > with m = 0) the term F],^^ dominates the term 
Fj^^O' which is only proportional to A° there. The application of this result to the late time 
tail will be discussed below. 

We now turn to the terms k >2, and show by mathematical induction that for all > 1 
there exists a solution admitting the form 

Fi^,{r)=A-^hk{r) (100) 

[with hk{r) being functions analytic at the EH], and thus satisfying the EH-regularity con- 
dition (160|) . This form was already verified in the case /c = 1, for which we also showed that 
/ii(r+)^0. 

Let us assume that Eq. ( |100| ) applies for all 1 < k < k' where fc' > 1 is arbitrary, and 
show that it is also valid for k = k'. Substituting the form ( 100| ) into Eq. ( ^3] ) we find 
Zl., = A~''^hk/{r) where hk'{r) is analytic at the EH. We then have 



V, X/ [0+(rO]2 [A(r')]^+i " ^"^^n {r)h,,{r) 



7fc'0+(r)ln^+ + A-^+i/ifc,(r), (101) 
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where hk'-, hk', and hk' are analytic at the EH, and 7^/ are constant coefficients. (Here we 
have taken r2 = t+.) For Fl, we thus obtain 

FUr) = at>^t + K' [^"'(1 + + ■■■)+ /3 In(^) 0+] + 7fc' In z+ 0+ + A-'+%ir). (102) 

This result is analogous to Eq. (|96|), only here we have not ruled out the possibility that 
some of the coefficients jk' may vanish (for k' > 1) . 

For the regularity condition to be met, we must now have b^, = —'Jk'/P- Consequently, 
Eq. ( [1001 ) is recovered for k = k' as well. By induction, then, we conclude that Eq. (|100|) 
applies to all > 1. 

Recall that it is possible for some of the coefficients 7^ (with A; > 2) to vanish. If, for 
a certain k, 7^ happens to vanish, then, to maintain EH-regularity, one must have 6^ = 0. 
In that case Eq. ( |102|) yields for Fl a divergence rate slower than A"'' at the EH (for this 
specific k). This, of course, does not contradict Eq. ( |100D , which should be regarded as 
merely setting an upper bound to the divergence rate of the functions -F'fc>2- It is only for 
F^^^ that we verified the actual asymptotic behavior oc A~*, by showing hi{r^) 7^ 0. This 
information, however, would be sufficient for the late time analysis at the EH. 

The above results, concerning the behavior of the physically-regular functions Fl{r) at 
the EH, are arranged in Table I. The table shows the leading-order forms of these functions 
for the various cases studied above. Below we use these results to discuss the late time 
behavior of the fields \E''' along the EH. 



the case 


Fk=o 


Fk=i 


Fk>2 


am = 0, s > 


A" 




< A"" 


all other cases 


A-" 6*""^^+^* for all A; > 



TABLE I. The asymptotic behavior of the "physically regular" functions ^^(r) at the EH. Presented 
are the leading-order in A forms of these functions, for the various cases studied in the text. For axially 
symmetric (m = 0) modes the asymptotic behavior depends on whether s > or s < 0, as discussed in the 
text. Note that these asymptotic forms arc in all cases independent of the multipolc number I of the modes 
under consideration. 



B. The late time tail along the EH 

When discussing the field behavior along the EH, it is most natural to refer to the 
"physical field" = A''\E'*, which, by construction, is a linear combination (with regular 
coefficients) of the regular Weyl or Maxwell components. Expressing the LTE in terms of 
and using Eq. (^8|), we find 

00 

fc=0 l,m 

Here, the factor in the square brackets is the function Fl{r) which, by the above construction, 
is regular at the EH (for all k). Recall also that the angular dependence here is EH- regular, 
and that the f-dependent factors take finite values at the EH. Thus, each of the terms in 
the sum over k in Eq. ( |103|) is indeed physically regular at the EH. 



29 



Now, at large the field should, potentially, be dominated by the k = term in Eq. 
( |103|) . For the s < fields and for m ^ modes of the s > fields, we find from Table 
I that at the EH itself the factor in the square brackets in Eq. ( |103D admits [ oc const 
for all k > 0. Therefore, in these cases, the late time decay of along the EH is indeed 
dominated by the k = term, with other terms smaller by factors of 1/v. Let us denote by 
yjysm part of \E'^ which includes all multipole modes I of a given m. To leading order in 
1/v, we then find for all modes m of the s < fields, and for non-axially symmetric (m 7^ 0) 
modes of the s > fields, 

^sm^^ = ^_^) = ^ ^^""{6, (^+) 6^"^^+^ v-^°. (104) 
I 

The situation is different in the case of axially symmetric (m = 0) modes of the s > 
fields. Here, the factor in the square brackets in Eq. ( |103| ) vanishes at the EH like A* for 
k = 0, whereas for = 1 it is finite. Hence, in this case, the k = 1 term dominates the A; = 
term. For each of the k > 2 terms the above factor is at most finite (and may even vanish 
for some k); hence these terms are negligible with respect to the k = 1 term at large v. We 
conclude that for axially symmetric modes of the s > fields the late time behavior along 
the EH is dominated by the k = 1 term: 

^">°'™=°(r = r+) = ^ 6+ Y'^'"'=°{e) v~'''>+\ (105) 

I 

to leading order in 1/f (here, the coefficient bi is also / dependent). 

At this stage, we still do not know the mode composition of the above leading-order 
tails; neither can we tell what the power index k^ and the amplitude coefficients are. These 
pieces of information will be obtained below, by matching the LTE to the form of the late 
time field at null infinity. 

However, one feature of the behavior along the EH is already manifested in Eq. ( |104| ) 
above: Non-axially symmetric (m 7^ 0) modes of the fields do not exhibit a strict power-law 
decay along the EH; rather, the amplitude of the power-law tail oscillates along the null 
generators of the EH, with an (advanced time) frequency mf2+. This phenomenon was first 
observed by Ori ||2^, and was further analyzed in Refs. pT| , p5| , pi . 



We comment that the above oscillations are not manifested when using the Kerr coor- 



dinate (f, defined by dip = dip + {a/ A)dr, instead of (cf. Ref. [|T8l, which adopts the 
coordinate (p>). Both coordinates are regular at the EH; however, the horizon's null gener- 
ators are lines of constant but varying (^.|^ Note that the oscillation of the scalar field 
along the horizon's null generators is a coordinate-independent phenomenon. 



To see that, we point out that, in the {v, r, 9, (p) system, the null generators of the EH are lines 
oi 9 = const, r = r+, = 2a A, and v = 2(r^ + o?)X, where A is an affine parameter along the 
generators — see §33.6 in Ref. [39|. Thus, along the null generators we find = O+u. Now, at the 
EH the two coordinates (p and are related by dip = + , from which we conclude that 
dipj^ = along the null generators of the EH. 
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VII. LATE TIME BEHAVIOR AT FIXED r > r+ 



In this section we obtain the global late time behavior of any of the modes of the field 
at any fixed value of r > r+. This task is to be accomplished in two steps: First, we 
consider the mode / = Iq (for each given s and m), which in Sec. |T| we found to be the 
(single) dominant mode at null infinity at late time. By evaluating the form of this mode of 
the LTE at null infinity, and comparing it to the form obtained independently in Sec. Ill, 
we derive the unknown power index ko, as well as the amplitude coefficient of the late time 
tail at fixed r. Provided with the value of ko, we can then carry on and, in the second step, 
obtain the tail form of all other modes at fixed r > r+. In particular, we then find that the 
single mode Iq dominates the behavior of the field \E''"™ also at fixed r. 



A. The mode / = / 







Substituting v = t + in Eq. (0), we find that along any r = const > r_|_ world line, 
at the late time limit t ^ \r^\, the behavior of the mode v]/'^'om, ^-^^^ each given s and m) is 
described by 

= a|)" y^'o™(^, (^) 0j,o(r) t-'^o (106) 

to leading order in \r^:\/t, where the second equality is due to Eq. (^Uj). To obtain the 
unknown power index k^ and the amplitude coefficient Oq , we now match the LTE to the 
form of the field at null infinity, as derived in Sec. |T^. 

In order for Eqs. (|48D and (0) to agree at null infinity for the mode I = Iq, we must have 

at scri+ 



.k=0 



Here, V'lte denotes the time-radial part of the mode Iq at null infinity, as calculated from 



the LTE, whereas the expression on the right-hand side (RHS) is the one derived in Sec. |T1 
using the iterative expansion scheme. As it turns out (see below), all terms k of the sum 
on the LHS of this equation contribute in the same order of magnitude at null infinity, and 
thus should be all summed up when evaluating '^/'lte. To that end, we first need to obtain 
the large-r asymptotic form of all functions Fj°{r). 

Starting with A; = 0, we have from Eq. (H), F^^""" = ao°0j?. Using Eqs. (|66|) and (|72D we 
find the large-r asymptotic form 

F^t^ir > r+) ^ r'o-^ (108) 

where the constant coefficient is given by 

^lo -Jo^i for s > with m = 0, , . 

7o - «o X I ^km ^^Yiei cases ^ ' 

[with the coefficients B'-/"^ given in Eq. (|6^) ]. 
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To analyze the functions i^fc>i we use Eq. (pO|), with the coefficients bj. taken to be the 
ones determined above by EH-regularity considerations (e.g. bk = for m 7^ 0). We now 
show by mathematical induction that for all k, these functions admit the large-r asymptotic 
form 



F^''(r »r+) 



(110) 



(to leading order in r), where 7^° are constant coefficients. To that end, we assume that Eq. 
( 11101) applies to all k < k' (where k' is an arbitrary integer greater than zero), and verify its 
validity to k = k'. 

First, we must calculate the large-r asymptotic form of the source function Zj^',. From 



Eq. 



we obtain, to leading order in r, 

Zj3 = 2{k, + F - 1) \r\dFiU/dr) + rF^?_; 



— Ik'^ ) 
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where = 2{ko + k' - l){lo - s + A;')7fc'-i are constants.^ 

Substituting this leading-order form in Eq. (|8^ ) and performing the double integration, 
we obtain, to leading order in r, (f)]^ = ^'■Oj.io-s+k ^ ^j^]-^ 



k'{2lo + k' + l) 



2{ko + k' -l){lo-s + k') 
k'{2lo + k' + l) 



Thus, for all k' > 0, the contribution from (j)}^ to PjS at large r [via Eq. (|80D] dominates the 
contribution from the homogeneous solutions, which is at most ~ r^o-s^ therefore have 
fI° = at large r; hence Eq. ( |110| ) is satisfied for k = k' as well. We also know by Eq. 
( |108|) that Eq. (|110D is valid for A; = 0. This, by mathematical induction, verifies Eq. (|110|) 
for all > 0, with the coefficients 7^" given by 



^0 

7fc 



(2/o + l)!7o° 
(A;o-l)!(/o-s)! 



2^(A;o + fc-l)!(/o-g + fc)! 
A;!(2/o + A; + l)! 



;il2) 



We are now in position to evaluate the sum over k on the LHS of Eq. (|107|) at null 
infinity. Substituting for Fk{r) [using Eq. ( 110 )], and recalling r = = {v — u)/2 at large 
r, we obtain 



LTE 



00 

lo-s-ko lo 



fc=0 



iQ—s+k 



;ii3) 



To evaluate the sum of this power series at null infinity, we make use of the auxiliary identity, 
valid for |g| < 1, 



Actually, the source Zj^, contains also contributions from other modes (/ = /q + li^o + 2). 
However, as we show later in this section, such contributions are negligible at large r, and do not 
affect the asymptotic form ( |111| ). 
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^ {ko + k - ly.ji - s + k)\ ^ 



k=0 



k\{2l + k + l)\ 



J-s \ 



(fco-2«-2) 



(114) 



where the derivatives on the RHS are with respect to q. [To prove this identity, insert 
the power expansion = q^~'^ J2T=o l'^ i^^o the RHS.] We now make the substitutions 
q — >• {1 — u/v) and / — > /q- At f m (recall that at null infinity v ^ oo whereas u 
takes finite values), the expression on the RHS of Eq. ( |114[ ) is then dominated by the 
term resulting from k^ — I — s — 2 differentiations of the factor (1 — q)~^ — which yields 
{kQ-lQ-s- 2)\{u/v)-'^^''-^^-'-^\ Substituting this result in Eq. (|rT3|), we obtain 



LTE 



2-'°+-(2/o + l)!(fco-/o-.-2)!7^° 



(fco-l)!(/o-s)! 
Comparison of Eqs. ( |1(J7| ) and ( |1 1 5| ) finally yields 

kf) = 21q + 3, 

and also, with the help of Eq. ( |109|) , 



u 



{ko-lo-s-l) 



, , 2'o-"+i(/o + 1) 



lo-s + 1 



l/Bbr~^ for s > with m = 0, 
llB^^"^ in all other cases. 



(115) 



(116) 



(117) 



The parameter /cq derived above is, in view of Eq. ( |106| ), the power index of the /q 
mode's tail at fixed r > r+ for each given s and m [recall /q = max(|s|, \rn\)\. The coefficient 
Oq describes the amplitude of this tail, with Eq. ( |1171) relating it to the amplitude of the 
/o-mode's tail at null infinity (which is ao°/o°). 



B. The modes I > Iq 



We now turn to analyze the behavior of the other modes, I > Iq, at fixed r. Here, the 
coupling between modes in the Kerr case shall appear to have a crucial effect on the form 
of the late time tail. To discuss this effect, it is most instructive to first consider a situation 
without coupling. Thus, at first, we shall "turn off" the interactions between modes by 
ignoring for a while all terms in [Eq. (^3])] which couples the mode I to other modes. 
(This will qualitatively describe the situation in the Schwarzschild case.) Then, in the second 
part of the following discussion, we restore the coupling (by taking into account all terms in 
Zl), and discuss its important effect on the late time tail of decay. 



1. The case with no coupling between modes 

Considering a mode / > Iq and ignoring its interaction with other modes, we can follow the 
same calculation as for the Iq mode, and obtain at null infinity iI>ite ^ Oq f 
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in a full analogy with Eq. ( |115|) .P^ However, since for any I > Iq we have ko — I — s — 1 = 
21q — 1 — s + 2 < l — s+2, this result cannot match the known u power at null infinity, 
unless flg = 0. The "boundary condition" at null infinity thus dictates the vanishing of Oq, 
and thus of Fq, for all / > Iq. This, of course, means that the modes / > Iq will decay faster 
than t~^° at fixed r. 

This argument demonstrates that in order to determine the form of the late time tail at 
fixed r, it is necessary to find out, for each given mode I > Iq, what is the smallest k value 
for which Fl does not vanish. Denoting this value by k{l), we find that the tail of any mode 
/ is given hj t'^"-''^^^ = t-2«o-3-fc(0_ 

What is k{l) then? We shall answer this question by matching the LTE at null infinity 
for each mode I > Iq. First, we must obtain the asymptotic form of the functions F^ at large 
r. By definition of k we have Fj^^j^ = 0. If no interaction occurs between various modes, 
then the function -F'r must be a static solution (the regular one): Fjr = al'rcj)^ = a'rr'~* at 
large r, where the constant coefficient does not vanish by definition of k. Following the 
same calculation as for the Iq mode (again, with the interaction between modes ignored for 
a while), we obtain for each k > k the large-r asymptotic form 

^k = 7ir'-^+^'"'\ (118) 

where 7^ are constant coefficients. Hence, for each mode I > Iq we find, to leading order in 

u/v, 

00 

4te ^Y.^kri-'-^'-'v-''-' = v-'-^'^^^'' X f{u/v), (119) 

k=k 

where / = E^fc7fc[(l - m/i;)/2]'"''+^"^ is a function of u/v only.0 

Now, we know independently from Sec. that any I mode decays at null infinity with 
a tail of the form oc ^-2s-iy~('-s+2)_ gy comparing the v power in Eq. (|119|) to this form, 
we find (setting /cq = 2/o + 3) that the function f{u/v) must admit the asymptotic form 
f{u/v) = [u/v)~'^^°^^^^^'^~'^. Then, by comparing the power of m, we finally obtain 

k{l) = 2{l - Iq) (no coupling; SBH case). (120) 

We conclude that for each mode I > Iq, the form of the field at null infinity (acting as a 
boundary condition) dictates the vanishing of all term of the LTE with k < 2{l — lo). The first 
nonvanishing term, the one with k = 2(/ — Zq), exhibits the late time tail t-^o-'i{i-io) _ ^-21-3 _ 
Other terms, the ones with k > 2{l — I) , decay faster at late time. Therefore, substituting 
V = t — in Eq. (|50D, we find for each of the modes / > Iq, 



Recall that in obtaining Eq. ( |115| ) for the Iq mode, the interaction between modes was not taken 
into account (this will be justified below). Thus, the analogy with this case is straightforward: one 
only needs to replace Iq with / in Eq. (|115D . 

Here we do not sum over k in an explicit manner, as we did for the Iq mode. Rather, we use a 
simpler argument, which is yet somewhat less rigorous. 
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"^'^""{t > |r,|) oc ¥'^""{9, if) (Plir) t'^'-^ (no coupling; SBH case) (121) 
to leading order in \r^\/t. 

2. The effect of coupling between modes 

Let us now "turn back on" the interactions between modes, and consider their effect on 
the late time tail. For each of the functions F^, the inhomogeneous part 0^^^ in Eq. (|80D 
now contains contributions not only from the functions -F^.i and Fl_2 but also from -F^^2' 
and (for s ^ 0) from F^^^ and -^^^2- ^'^^ example, the function Fj^!^^^ admits (for s ^ 0) 
a nonvanishing source oc aFjj'^^ = aoQ^^". Since we have Og 7^ (by definition of /cq), 
then, necessarily, the function F^Pj]^ is nonvanishing. Thus, for the mode I = Iq + 1 we find 
k = 1, implying a late time tail of the form t~^o~'= = t~2«o-4 'jj-^jg jg different than in the 
"no-coupling" situation, in which for the I = lo + 1 mode we had k = 2 [Eq. (|120|) ], leading 
to the tail t-2«o-5 j-gq^ ( |121| )]. We may summarize the result in this example by saying that 



the interaction "excites" the mode + 1 already at = 1, whereas the boundary condition 
at null infinity "excites" this mode only at = 2. We arrive at the conclusion that for this 
mode the effect of interaction dominates the late time behavior! 

The generalization of this result to all I > Iq is straightforward: For each I > Iq, the 
function Fj^ [recall that k{l) is, for a given mode /, the smallest k for which Fj^ is nonvanishing] 
serves as a source, via Z^, to -F/"*"^ and (for s 7^ 0) to Fl~^] and F^^i. These three functions 

' k' k+2 ^ ' ' fe+1 fe+2 

are then necessarily nonvanishing. Since we have kil = Iq) = 0, we obtain in the case s 7^ 0, 
k{l) = l — lo for all / > lo- In the scalar field case (s = 0) the interaction couples only between 
next-to-nearest modes. In this case, the mode I = Iq + 1 (for example) is not excited by 
interaction, but rather by the boundary condition at null infinity, yielding k{l = Iq + 1) = 2 
[see Eq. (|120|) ]. In general, for s = we thus find k = I — Iq for even / — Iq, and k = I — Iq + I 
for odd / — Iq. We can express the results in all the above cases by writing 

k(l) =1 -l^-^q (with couphng), (122) 

where 

_ J 1 for s = with odd / — Iq, , . 

^ 1 otherwise. 

By comparing Eqs. ( |120|) and ( |122|) we find that interaction excites any of the modes 
I > Iq at smaller k than do the boundary conditions at null infinity. In other words, it is the 
interaction (rather than the boundary conditions at null infinity) which first excites any of 
these modes, and thus determines the form of the leading order tail at late time. This tail 
shall admit the form t~^°~^ = t-'-'o-3-q_ 

More precisely, for each given s and m, and for each of the modes I > Iq, we find from 
Eq. (pOD, to leading order in |r^,|/t, 

^''""{t > |n|) oc F^'"(e, ip) FLj_,„+g(r) (Kerr case). (124) 

Note that for each / > Iq (for the scalar field — for each I > Iq + 1) the late time tail decays 
s/oiyer than the corresponding tail of the same mode in the Schwarzschild case, Eq. ([1211) . 
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These slowly decaying tails are produced by interaction, as discussed above. On the other 
hand, the mode Iq (for each s and m) , whose form is not shaped by interaction, exhibits the 
same decay rate in both the SBH and the KBH cases: 

There now remains an important question to deal with: In what way does the matching 
of the LTE at null infinity change under the effect of interaction between modes? Does the 
LTE remain consistent with the boundary condition there? In particular, we would like to 
show that the calculation made above for the Iq mode is still valid even when the interaction 
is taken into account. 

To answer these questions we must first examine the large-r asymptotic form of the 
various functions F^, this time taking into account also the effect of interaction between 
modes. As we recall from Eq. (|80D , each function Fl contains a "homogeneous" part, cii0^ + 



b^k4>\rj "inhomogeneous" part, [0y*^. The homogeneous part vanishes identically for 

all k < 2{l — Iq) [by virtue of Eq. ( |12CI|) ], and behaves as r'~* at large r for k > 2{l — Iq). The 
inhomogeneous part vanishes identically for all k < I — Iq + q [by virtue of Eq. ( p.22| )]. Its 
large-r form can be calculated for each k and / using Eqs. (^) and (p4D, given the large-r 
form of all functions Fl, which serve as sources to Fl 



k- 

Using Eqs. (p3D and ([8^), one may formulate two practical calculation rules: 



(i) If the function admits the large-r form ~ r^ (where p is some positive power 
index, and a '~' symbol represents the asymptotic form to leading order in r), then the 
contribution to Fl due to the term in involving Fl_^ would be of order ~ t-p+i. For 
example, the source function -F^Lq (which is an EH-regular homogeneous solution and thus 
admits ~ r'o-s^ induces on -F^Li (through Z^^l) a nonvanishing contribution of order ~ 

(ii) The contribution to a function Fl due to interaction with a source function -F^/"^' 
admits the same large-r form as the source function itself. This is the situation with all 
interaction sources of Fl and, also, with the source function Fl_2. For example, the source 
Fjl^Q ~ r'""** induces on -F^^^ (through Z^~2^^) a nonvanishing contribution of order ~ r'°~**. 
It also, for example, contributes in order ~ 7-'°"^* to ^^^2 (through 2'[~2)- 

With the above two "rules of thumb" one can now inductively construct expressions for 
the large-r form of each of the functions Fl, starting with F'feLo ~ r'°~'*. This easily yields 
Pl ^ ^{io-s)+k-(i-io) ^ namely, 

Fl{r > r+) ~ r2'o-^-^+*^ (125) 

for each / > Zq and k > k{l) = l—lo+q^ This is the same asymptotic form as obtained for the 
functions Fl in the absence of interaction, for k > 2{1 — Iq) — see Eq. ( |11^ ) with k = 2{1 — Iq). 
[The important difference is that without interaction the terms I — Iq + q < k < 2{l — Iq) all 
vanish identically. 



In fact, this result, Eq. ( |125D , is not completely accurate, as the functions F^^^ turn out to 
involve logarithmic factors which complicate the situation. This logarithmic dependence will be 
discussed at the end of this section; meanwhile we shall ignore it to make the discussion more 
clear, and refer only to the power-law dependence of the functions Fl (which is not affected by the 
presence of the logarithmic factors). 
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We conclude the following: 

(I) For any mode /, the "inhomogeneous" part (~ j^'^io-i-s+k^ pi^ dominates its "homo- 
geneous part" (~ r'~'*) at large r for all k except k = 2(/ — /q), where both parts contribute to 
order ~ r'~* [recall that the homogeneous contribution vanishes identically for k < 2{l — lo)]. 
Thus, for each mode /, only one arbitrary parameter [the one belonging to the homogeneous 
solution at = 2{l — Iq)] is involved in the leading order form at large r — as we may expect. 
(For each I, this parameter is to be determined, in principle, by matching at null infinity.) 
Note also that the first non- vanishing function of each I (namely is always proportional 
to the parameter Oq , which originates from the mode I = Iq aX k = 0, and "propagates" 
through interaction to higher modes. [This parameter was determined above, Eq. ( |117| ), by 
matching the mode /q at null infinity.] 

(II) Whether the interactions are taken into account or not, one obtains the large-r form 
Fl oc r'^^o-i-s+k j^Qj, k > 2{l — Iq)], though with different proportion coefficients. That 
difference only affects the amplitude of the function '?/'lte [^-g- in Eq. ( |119| )], yet the matching 
of the modes I > Iq a.t null infinity, discussed above, remains qualitatively the same. 

(III) As to the effect of the interaction on the mode Iq: By Eq. ( |125D we have for all k, 
at the leading order in r, fI° ~ j.io-s+k^ This function has, in general, interaction-induced 
contributions coming from -F^^^, FI~^\, and -^^^2' which by Eq. ([L25|) and the above rule 



(ii) are of order ~ ^ ^^^^ ^ ^lo-s+k-A^ respectively. Thus, the large-r 

leading order behavior of the functions fI° is not affected by interaction with other modes. 
This result is valid only for the Iq mode (and, for s = 0, also for Iq + 1) which only interacts 
with modes of larger /. It justifies ignoring the interactions when evaluating the behavior 
of the mode Iq at null infinity, as we did above. Particularly, the values derived above for 
the parameter fco [Eq- (|116|) ] and for the coefficient ag [Eq. (|117|) ] remain valid also with 
interaction between modes taken into account. 

Finally, it should be commented on that in the above discussion (regarding the / > 
Iq modes) we have ignored a certain complication for the sake of clarity: Actually, the 
integration in Eq. (Q) produces, for certain values of k and /, a logarithmic dependence 
on r. This leads to an asymptotic form of Fl which is not strictly a power law [as in Eq. 
( |125|) ], but, in fact, having the form Fl ~ j,2io-i-s+k ^ (Inr)'^. It can be shown that for 
each / and k, the logarithmic power L may take only integer values between and / — Iq. 
Particularly, we find no logarithmic dependence in all functions F^ belonging to the mode 
Iq, and thus no modification is required in the above analysis for this mode. Also, we can 
show that L = for all modes F| (i.e. for the first nonvanishing function F^, of each mode 
/), and therefore the functions Fl in Eq. ( 124 ) exhibit no logarithmic dependence at large r. 



Finally, the matching of the / > Iq modes at null infinity, discussed above, was based merely 
on the power-law dependence of the functions Fl{r), which is not affected by the presence 
of the logarithmic factors. However, it is not clear to us whether the logarithmic factors 
themselves properly match at null infinity [recall that in Sec. |T| the logarithmic dependence 
(in u) of the modes I > Iq at null infinity has not been fully investigated]. This question 
remains open. 
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VIII. SUMMARY AND DISCUSSION 



In this paper we have explored analytically the late time decay of the Newman-Penrose 
scalars \E'* (representing scalar, electromagnetic, and gravitational perturbations) in the 
background of a realistic Kerr black hole. Our analytic method can be summarized as 
follows: We assume that at late time each of the fields \E'* admits the late time expansion, Eq. 
(^). This reduces the master perturbation equation to a hierarchy of ordinary differential 
equations for the radial functions Fl{r). The homogeneous part of each of these equations 
is just the static field equation in Kerr spacetime, to which there exists an analytic basis 
of exact solutions. In addition, for each / and A; > 0, each of these equations possesses 
an inhomogeneous part depending on functions Fki<:k (including functions which belong to 
other modes). Using the Wronskian method we can then explore, in an inductive manner, 
the general solution for each of the functions F^. Each of these solutions contains, in advance, 
two unknown parameters. One of these parameters is determined by regularity requirements 
at the EH. The other parameter is determined by the form of the field at null infinity (serving 
as a boundary condition). To obtain the behavior of the fields at null infinity, we apply the 
iterative scheme, as described in Sec. [m| . 

Following is a summary of our main results. These results are valid in the most realistic 
initial setup of a compact pulse composed of all multipole modes (and in particular, the 
lowest radiatable mode Iq for each value of m; below we also briefly discuss the more special 
case where this mode is missing). 



Tail form at fixed r > r+ 



Along any worldline of fixed r outside the KBH, each specific mode /,m decays at late 
time with the tail 

^''"'(t > |n|) ocr('+'°+=^+'?) for each Z>/o = max(|s|,|m|), (126) 



to leading order in \r^\/t [see Eq. ( |124| )]. Recall that g = 0, except for s = with odd / — Iq 
in which case q = 1- 

The most dominant modes of the overa// field are those with / = |s| and — |s| < m < 
\s\. From Eq. ( |106| ) we find, to leading order in |r*|/t, 

l«l 

^^(t>|r,|)= ao"'''>"'''=''''™(^,<^)0t"'''(r)r(^l"l+^) (overall behavior). (127) 

r?i=— |s| 

Here, the function 0j.(r) is the physically regular static solution, whose exact analytic form 
is given in Eqs. ( |62a| ), (^^, and ( |79a| ) , corresponding, respectively, to the case s 7^ with 
m 7^ 0, the case s 7^ with m = 0, and the case s = 0. The constant coefficient oq (which is 
also m dependent) is related in Eq. ( |117| ) to the amplitude of the leading-order tail at null 
infinity, which, in turn, is expressed as a functional of the initial data function — see Eq. (^5]) 
and the Appendix. Note that Eq. ( |127D constitutes an exact analytic expression (accurate 
to leading order in \r^\/t) for the late time behavior of the fields \Ef*, valid anywhere at fixed 
r > r+. 
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The power law indices predicted in Eqs. ( p.26| ) and (|127|) agree with those obtained by 



Hod 1 24 1 at fixed r ^ M (the result by Hod refers only to this asymptotic domain). The 
result in Eq. ( |127| ) has support from numerical simulations (in 2+1 dimensions) by Krivan 
et al. — see Ref. for s = and Ref. |19| for s = —2. Also, the validity of our prediction, 



Eq. ( |126| ), has recently demonstrated numerically by Krivan 



Tail form along the EH 

It is most natural to express the results at the EH in terms of the "physical" fields 
^ A^\E''', which are related through an EH-regular transformation to the components of 
the Maxwell and Weyl tensors (see the discussion in Sec. |IV| ). By virtue of Eqs. ( p.04| ), (|105D , 
and ( |122| ), we find each specific mode /, m to decay along the EH with the tail 

jsim( I \^ j y-{i+io+4+Q) fors>0, m = 0, , . 

W (t; > |r,| j oc I ^-ii+io+3+q) in g^ii other cases, ^ ^ 

(to leading order in \r^\/v), with defined in Eq. (0). Note that for the s > fields 
the axially symmetric (m = 0) mode decays faster than other modes. Consequently, the 
late time behavior of the overall field is dominated by the non-spherically symmetric, 
m 7^ 0, modes. These modes oscillate along the null generators of the EH with (advanced 
time) frequencies mfi+. We find, to leading order in \r^\/v, 

Ejl|=i 4='^' l^^''=l^l'" (0, ^+)e^™^+^ v'(W)^ overall s > field, 
^ E 1=0 4^''' overall s < field, ^ ^ 



where the regularized azimuthal coordinate (p+ is the one defined in Eq. (p5|). Thus, the late 
time behavior of the field vl/^^" along the EH is characteristically oscillatory. On the other 
hand, the behavior of the scalar field {s = 0) is characteristically nan- oscillatory, whereas 
the field involves both oscillatory and non-oscillatory modes. 

As recently discussed by Ori ||2^, the characteristics of the late time decay along the 
EH — both the value of the power index and the oscillatory nature of the waves — have 
important implications to the structure of the infinite blue-shift singularity at the inner 
horizon of the KBH. This singularity is related with the behavior of the ingoing component 
\l/o = of the Weyl perturbation. As it turns out [^, this singularity is generically 

oscillatory. 



Tail form at null infinity 

For each specific mode I, m, the analysis in Sec. ^ predicts at null infinity a late time 
tail of the form 

^"'™(U > M) OC r^''"(^, ^) V-^'-^ ^ ^^3Q) 

to leading order in Uq/u and in M/uq. (We assume here that the initial pulse is emitted 
at large distance, so that —uq ^ M.) In this expression, "[l.d.]" represents a possible 
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logarithmic dependence of the form ln^[u/(r+ — r_)], where the power L is some positive 
integer. Such a logarithmic factor does not occur for the dominant modes Iq of each m, and 
also for all modes of a scalar field (s = 0). Our analysis indicates that logarithmic factors 
do occur for the less dominant modes (/ > lo) of s ^ field; however, this point was not 
studied by us in full detail. In any case, for each given s and m, the dominant late time 
decay at null infinity is described by Eq. (^Hf), in which no logarithmic factors occur. 

The overall field \E''* is dominated at null infinity by the modes with / = |s| and — |s| < 
m < |s|: 

^'(m>M)= ao~^'^ Io~^'^Y''^=\'\''^{e,ip)v~^"'^u-^\'^"+^^ (overall behavior) (131) 

m=—\s\ 

to leading order in uo/u and in M/uq. Here, /q"'*' is a functional whose construction from 
the initial data is described by Eqs. ( ^5]) and (|36|). 

The power-law indices given in Eqs. ( |130|) and (|131|) are in agreement with Hod's results 
23|,|2^, though Hod indicates no logarithmic dependence for any on the modes. 



Non-generic initial data 

We now briefly discuss the case where the initial pulse is of a non-generic mode compo- 
sition, such that (for given s and m) it does not contain the mode / = Iq. For example, what 
can we say about a case in which, for a s = ±2 field, the angular dependence of the initial 
pulse is that of a pure mode m = 0, Z = 4? 



The calculation scheme presented in Sec. Vll, based on the LTE, allows one to obtain 



the power index of the tail at fixed r regardless of the initial setup, provided only that the 
power index at null infinity is known. Suppose that, for a specific initial setup and for a 
certain m, the most dominant mode, Iq = max(|s|, \m\), falls off at null infinity with a tail of 
the form "^"^"^ oc u~^. Then, from Eq. (|115|) (whose derivation does not involve any reference 
to the details of the initial data) we must have w = kQ — Iq ~ s — 1, where is the power 
index of the tail at fixed r. Therefore, for any initial mode composition, there exists a simple 
relation between the power law indices at null infinity and at fixed r. Symbolically, we may 
write (for each m) 

vl/^-^cxtz- ^™t>.+ cxt-(-+'«+^+^), (132) 

which is valid for any initial mode composition. 

The main challenge, then, is to obtain w, the power index at null infinity. Our iterative 
scheme, presented in Sec. |II|, provides a formal way for accomplishing this task; however, in 
the case of a non-generic initial mode composition this technique becomes less practical, for 
the following reason: Consider, for example, an initial pulse of scalar radiation, composed of 
only the mode m = 0, Z = 4. Then, the function tpn=i (namely, the first-order iteration term; 



we use here the notation of Sec. p| ) would contain only the three modes / = 2, 4, 6. The 
most dominant mode of the overall field, that is Z = 0, would be excited only at n = 2 (and, 
as we suspect, will gain its typical power-law form only at = 3). Thus, in this example, it 
would require one to go through at least three successive iteration stages in order to recover 
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the tail form of the dominant mode. (Recall that in this paper we only discussed the first 
iteration; the second iteration already becomes very complicated for analytic treatment.) In 
general, as larger is the difference between the initial mode and Iq, as greater becomes the 
number of iterations required to extract the exact tail form of the dominant mode. It is only 
for the generic case discussed in this paper that a single iteration suffices for this goal. 

The case of non-generic initial data (specifically, the case of any pure initial mode) was 
studied by Hod in Refs. and However, recent numerical experiments by Krivan pB[ 
show disagreement with Hod's results in this case. Further work is needed, both analytic 
and numerical, to clarify this point. 



Final remarks 

We recall that in this paper we have considered only non-extremal, |a| < M, Kerr BHs. 
Clearly, the extremal case needs to be analyzed separately [note, for example, that Eqs. 
(p!9|)-(|2T|) cease to be valid in the case \a\ = M]. A basic property of the effective potential 
in both the SBH and the non-extremal KBH spacetimes — its exponential decay towards 
the EH (with respect to r*) — is no longer valid in the extremal Kerr case: rather than 
V{r) oc e'^'^+^* for \a\ < M near the EH, one finds V{r) oc r^^ for \a\ = M [where in both 
cases the tortoise variable r^, is defined through the differential relation dr^/ dr = {a^+r"^) / A] . 
Consequently, basic parts of the analysis presented in this paper may fail to apply in the 
extremal case. In particular, the crucial assumption made in Sec. |T1|, that the late time 
tail at null infinity is exclusively dominated by waves scattered at very large distances, need 
not necessarily hold in this case: Here, strong contribution to the tail may occur also from 
back-scattering at small distances! To clarify the situation in the extremal case, a separate 
detailed analysis is thus required. 

Finally, we should comment on the limited significance of "multipole modes" in the Kerr 
spacetime: The spin-weighted spherical harmonic functions y^^'™ are not related here to 
an underlying symmetry group, as they are in spherically symmetric backgrounds. Conse- 
quently, the "multipoles" associated with these functions have no invariant meaning, but 
are rather related to (and defined through) the specific choice of the coordinates. Yet, the 
functions y*'"^(^, 99) with 6,(p being the Boyer-Lindquist coordinates are signified as the 
natural basis for our purpose, because, at the late time limit, the field equation becomes 
separable in terms of these functions. (This separability is manifested in the k = term of 
the LTE, which exhibits no coupling between the various modes l,m.) Note also that the 
late time behavior of the overall field is in all cases governed by a pure mode I. 
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APPENDIX A: CALCULATION OF VLi AT NULL INFINITY 



In this appendix we calculate the "Ist-order" iteration term ip'n^i at null infinity, at 
large retarded time (that is, for v oo with fixed u ^ M) . In the case s > we shall give 
full details of the calculation. For brevity, the case s < (which turns out slightly more 
complicated to analyze for technical reasons — see below) will be discussed in less detail. 

The starting point for the calculation would be Eq. (|43|), in which the source function 
Sl^^i is given by Eq. (|31| ) as a function of various modes of 4'n=o- The radial functions SV{r), 
6R{r), and K{r) involved in the expression for 5"^=! are implicit functions of = (f — u)/2 
(because r is an implicit function of r*). We can expand each of these radial functions in 
powers of 1/r*. By virtue of Eqs. (^) and (^8]) we then obtain the leading order forms 



5V{r) 



M - imsa + l{l + 1)M [2 In 



1 



2rl 



O 



(Inr*) 



(Ala) 



and 



5R{r) 



m4-3M + 2Mln(- 



2r 



O 



(Inr*) 



(Alb) 



with the asymptotic form of K{r) given in Eq. ( p5c| ). 

The various terms in the source Sl^^^ contribute additively to V'i=i [via Eq. (^31)]. The 
analysis below implies that the dominant contribution to V'i=i at null infinity at large u 
comes only from the leading-order form (in 1/r*) of each of these source terms: Roughly 
speaking, each additional 1/r* factor in the source leads to an additional factor of 1/u in 
the contribution to ipn=i at null infinity. Hence, to calculate ipn^i to leading order in \/u, 
one may replace the actual functions 5V , 6R, and K{r), with their above asymptotic forms. 



Let us now consider the contribution to due to a source term of the form 



Sl^^{u',v';l'Pd) = (2n)-^ [a + /31n(2n)] 



(A2) 



where a tilde symbol over a quantity shall represent the ratio of that quantity and (r+ — r_) 
[so that = r^/{r^ — r_)]. This is the general form (to leading order in 1/r*) of all 
contributing terms in 5*^=1, with the integer numbers P, I', and d admitting the possible 
values P = 2, 3, = /, / ± 1, ±2 and d = 0,1, 2, and where a and /3 are constant coefficients. 

Let us denote the contribution of Sl^^^{u' ,v']l' Pd) to ipl^. 
Then, from Eqs. (^3]) and (|3^) we find 



at null infinity by i^n=i- 



X [a + I3\n{v' -u')]. 



u 



(n 



u 



' _ u'Y+i'-j+P 



(A3) 



Note that for s > 0, all derivatives with respect to u can be "taken out" of the v' integration, 



due to the factor (v' 



u 



d+s 



appearing in the integrand. This manipulation (which is not 



possible for s < 0) much simplifies the calculation in the s > case. For brevity, we 
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therefore continue from this point by concentrating on the case s > 0. Our calculations 
for s < 0, whose details shall not be presented here, yield the same qualitative results as 
those obtained below for s > 0, yet they are slightly more tedious (as they involve more 
complicated combinatorial expressions). Unfortunately, we could not figure out a way for 
treating the s < case in a simple manner as the s > case, though we think this should 
be possible. 

After "taking the u derivatives out" of the integration over v', we can now easily extract 
the u dependence of ipn=i' by transforming in Eq. ([A3| ) to the new integration variable, 
x{v') = {u — u')/{v' — u'). This yields 



(^ 



f du' [giW)]'^^^''^ X 

j=0 "'"0 



u — u 



i\i-i'-p+i+j 



aj {a + (3\n{u — u')) — f3jf3 | , 



(A4) 



where a,- and are constant coefficients given by 



{r -s + p-2-j)\{i + s) 



Jo 



{l + l' + P-l-j)\ 



(A5) 



[The coefficient aj is just the standard beta function, B{1' — s + P — 1 — j,l + s + 1).] 

Next, in Eq. ([A4| ) we integrate by parts j + d successive times with respect to u'. All 
resulting surface terms vanish due to the compactness of (m'), and one is left with 



I'-s 



^l-s+j+d 



[U — U 



i\l-V-P+l+j 



ocj {a + (3 ln(u - u')) - } 



(A6) 



Here, we have used the fact that the u' derivatives operate on functions of {u — u') only, 
to make the replacement du' — > —du- Evaluated at late retarded time, u ^ —uq, the last 
expression takes the form (accurate to leading order in u/uq) 



'o 



I'—s jl—s+j+d 

^ i rto,l-S+j + d 



+3 



aj {a 



+ /31nu) -I3jl3\] 



(A7) 



where /q is the functional constructed from the function 5'o'(u) according to Eq. (^). 
Finally, performing the multiple differentiation in Eq. ( |A7|) , we find 



jM'Pd jV 
Vn=l — -'O 



A 



WPd 



a + /? In m) + I3ri' 



WPd 



XV u 



s „-{l'-s+P-l+d) 



(A8) 



in which A^' ™ and rj'''''^'^ are constant coefficients. Here, the term proportional to A'' '^'^ 
contains only contributions which arise from all I — s + j + d derivatives in Eq. ([A7| ) acting 
on the power with non acting on In-u. The order of differentiation, I — s + j + d, 
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is in all relevant cases greater than the power index I — I' + j — P + 1. Hence, the only 
contributions to the coefficient A'' ™ arise when the power I — I'+j — P + 1 is negative, i.e. 
from the terms with j<l' — l + P — 2. Since the index j takes no negative values, we find 
that there would be no contribution to A'' ™ unless P >/ — /' + 2. Namely, 

For P >/ — /' + 2, the coefficient A"'^'' is given by 

,P>i-V+2 il + ^)'- 4r aVs( 1 ^i'-s+,+d il' -s + P-2-3)\{l'-s + P + d-2)\ 

^" -(I^l^o ^'^^^ (/ + /' + P-l-j)!(/'-/ + P-2-j)!' ^ ^ 

where j = min(/' — s, /' — / + P — 2) (in accordance with the above discussion) . An expression 
for the coefficient 77" ™ can also be written down explicitly, using Eq. ( |A7| ). It can be verified 
that this coefficient is non- vanishing for all relevant values of l,l',P, d. 

Equation (|A8D suggests that, potentially, one may expect logarithmic dependence to 



occur at the leading order tail at null infinity. In the following we show, however, that due 
to vanishing of the coefficient A'' ^'^ in certain cases, this logarithmic dependence is avoided 
as far as the most dominant mode of the overall field is concerned: this dominant mode shall 
appear to die off with a pure power-law tail. 

Note also that Eq. ( |A8| ) confirms our above assertion, that to leading order in 1/u, the 
contribution from each given term in S^^i (with given /' and d) comes merely from the 
leading order in l/r^: For higher-order terms in the expansion of the source there 

correspond larger values of P, leading to a faster decay in Eq. (|A8|) .H 

Using Eq. ( |A8| ) we can now analyze the contribution to V^^^i at null infinity at late 
retarded time, belonging to each of the various source modes /': 

a. Contribution of the mode I' = 1 — 2: For this source mode we have d = 2, and, at 
the leading order in P = 2 and /? = [see Eqs. (|3lD and (|25c| )]. Since for this order 
P </ — /' + 2 = 4, the corresponding A coefficient vanishes in Eq. (|A^). Thus, we find no 
contribution at all to "^^^i from the order O [{r'^)^"^] of this source mode. Turning next to the 
following order, with P = 3, we find again that A vanishes (as P < 4). However, at this order 

the logarithmic coefficient (3 of the source does not vanish [we have (3 = SMa'^C': — see Eqs. 

(pip and (25c)], and from Eq. (^) we find the nonvanishing, nonlogarithmic contribution 



V'i-ir' = ma^ci_ vp=3:L2 4"' (All) 

where we adopt the notation ^^=1 represent the late time contribution to ^n=i at null 
infinity due to the mode /'. Here, the symbol ' = ' stands for "leading order in M/u and in 

Uq/u". 



Higher-order terms in the l/r'^ expansion of the source would exhibit higher Inr'^ powers, 
leading to higher logarithmic powers in Eq. ( |A8| ); however, the u~^^ -s+P-i+d) pQ^g^ jg^^ would 
remain the same for the higher order terms as well, with P denoting the power of for each 
term. 
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b. Contribution of the mode I' = l — l: For this source mode there are two contributions: 
one with d = 1 [the one proportional to c'„ — see Eq. (0)] and one with d = 2 (proportional 
to C[_). Both contributions have P = 2 and /? = at the leading order in 1/r'^. Clearly, in 
view of Eq. ( [A8|) , the term with d = 1 dominates the contribution from this mode. Now, for 
the leading order, O [(r^)^^], we have P </ — /' + 2 = 3; thus the corresponding A coefficient 
vanishes. Since for this order we also have (3 = 0, one finds no contribution from this order to 

The dominant contribution would come from the next order (with P = 3), for which 
both coefficients A and j3 ■ r] are nonvanishing. Therefore, this contribution will contain a 
logarithmic dependence: 

^ 8iasM ci J^-i [A^='3;]=i (1/2 + lnu) + r]p=3j=i] u~^^-'+^\ (A12) 

c. Contribution of the mode V = 1: There are three terms in the source 5*^=1 which are 
not due to interaction with other modes [see Eq. (|31|)]. These are (i) the term proportional 
to 5V , for which P = 3 and c? = 0; (ii) the term proportional to 5R, for which P = 2 and 
d = 1] and (iii) the term proportional to Cg, with P = 2 and d = 2. Clearly, the 6V term 
and the 6R term contribute to ipn=i at the same order of whereas the contribution from 
the third term in smaller by a factor of 1/u. We thus concentrate on the first two terms, 
both of which have P>1 — I' + 2 = 2 already at the leading order in Hence, the 
dominant contribution to '?/'i=i would come from this leading order. From Eq. ([X^), using 
Eqs. (|Ala|) and ( |Alb|) , we obtain for the contributions of these two terms. 



^i^Udue to 5V) ^ -AM ll^v"u-^^-'+^^ 

X {Xp%^d=o [1 - imsa/M - l{l + 1) + 2l{l + 1) \nu] + 2l{l + l)r/^='3^rf^o} , 

(Al3a) 



<-^l(due to 6R) 



X 



-2Msll^v-'u-^^-'+^^ 
\i'=i 



,1 (ic[a/M - 3 + 21nM) + 2r/^ 



i-i 

P=2,d=l 



Now, from Eq. ( |A10| ) we find 



\p=3,d=o 



A 



i'=i 

P=2,d=l 



s{-iy-%i - s + + s) 

2/(/ + l)(2/ + l) 



21 + 1 



(A13b) 



(A14) 



Substituting these values in Eqs. ( [A13a| ) and ( |A13b| ), and adding up these two equations to 
construct the overall contribution from the mode /' = /, we find that the two logarithmic 
terms exactly cancel each other. (In the scalar field case, s = 0, each of the two logarithmic 
terms vanishes independently.) The leading order over- all contribution from the mode I' = I 
will therefore be nonlogarithmic, exhibiting a strict power-law tail of the form 



(A15) 
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d. Contribution of the mode I' = I + 1: For this mode, there contribute two terms in 
S^^^i, the ones proportional to c^,_ and to C^. These two terms have, respectively, d = 1 
and d = 2, and both have (at the leading order in P = 2 and (3 = 0. However, both 
coefficients X''j^2,d=i ^'^1^24=2 to vanish, resulting in the dominant contribution 
to ipn=i coming from P = 3. The contribution from the term proportional to c''_^ (which is 
oc n"^'"*"*"^) ln£t) dominates the one from the term proportional to C\_ (oc u^^'-~'^+^^ Intt), and 
one finds, in summary, 

^ g .^^^ ^i+i J^«^=i+i^^ (1/2 + In u) + rf'p=J^l^^ u~^^~'+'^\ (A16) 

e. Contribution of the mode I' = I + 2: This mode has d = 2 and (at the leading order 
in P = 2 and /3 = 0. However, the coefficient \^^2d=2-i well as Ap'z'g^J^g) turns 
out to vanish, resulting in the leading order contribution from this mode to ipn=i coming at 
P = 3 from the term in Eq. (|A8|) proportional to (3 ■ rj. Hence, this contribution would be 
nonlogarithmic: 

^i+2-.i ^ ^Ma'C'^^ Vp=3%2 (A17) 



Equations ( |A11|) , ( |A12|) , ( |A15|) , ( |A16| ), and (|A17| ) describe the various contributions to 
the tail of ipn=i at null infinity from all various source modes. These results are summarized 



in Eq. ( ^4t ) (in Sec. [HI D|) using a different notation for the amplitude coefficients. We point 



out that Eq. ([A8| ) , as well as all power- law formulas derived in this appendix, is also valid 
in the case s < 0, though with different amplitude coefficients. (As we mentioned above, 
we found these coefficients to be more complicated to calculate for s < 0; still, the same 
coefficients A found to vanish for s > 0, also appear to vanish in the case s < 0, which finally 
leads to the same power-law contributions.) 

The following are the main conclusions that can be drawn from the analysis in this 
appendix: 

(i) In general, the dominant contribution to the late time tail of a mode / of ipn=i at null 
infinity is due to the source modes /, 1 — 2, and (for s 7^ 0) / — 1. These contributions all have 
the form oc m^"^"*"^ (multiplied by a logarithmic factor in the I' = I — 1 case). Contributions 
due to the modes / + 1 and / -|- 2 are negligible. 

(ii) Consequently, for given s and m, the most dominant mode of the field ipf^i at null 
infinity is the lowest radiatable one, namely the multipole I = Iq = min(|s|, \m\). 

(iii) This mode (/q) admits no contributions from lower, I < Iq, multipoles, and thus, to 
leading order in 1/u, it is not affected by interaction with other modes. Equation ( |A15|) 



then implies that this mode (and thus also the overall field V'n=i) admits the strict late time 
power-law tail oc ^i-('o-s+2)^ with no logarithmic dependence. 
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